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(D ■ ABSTRACT 

Evidence favouring a Gaussian initial mass function for systems of old globular clusters 
' has accumulated over recent years. We show that an approximately Gaussian mass 

\ function is naturally generated from a power-law mass distribution of protoglobular 

T-H ■ clouds by expulsion from the protocluster of star forming gas due to supernova activity, 

^ ' provided that the power-law mass distribution shows a lower-mass limit. As a result of 

OO \ gas loss, the gravitational potential of the protocluster gets weaker and only a fraction 

■ of the newly formed stars is retained. The mass fraction of bound stars ranges from 

zero to unity, depending on the local star formation efficiency e. Assuming that e 

, is independent of the protoglobular cloud mass, we investigate how such variations 

affect the mapping of a protoglobular cloud mass function to the resulting globular 
cluster initial mass function. A truncated power-law cloud mass spectrum generates 
bell-shaped cluster initial mass functions, with a turnover location mostly sensitive 
1-^5 ' to the lower limit of the cloud mass range. Assuming instantaneous gas removal and 

a slope a ~ —1.7 for the cloud mass spectrum, we evolve the derived cluster initial 
Q i mass functions up to an age of 13 Gyr in a potential like that of the Milky Way. We 

( ' obtain a good match to the Old Halo cluster mass function, with a present-day mass 

\ mass fraction of clusters in the halo of 2 %, as is observed, with miow ~ 6 x 10^ Mq, 

■ niup > 5 X 10^ Mq, 5 ~ —2.9 and Tc — 0.025, respectively the lower and upper limits 
^ \ of the cloud mass range, the slope and the core of the power-law spectrum for the star 

• i-H , formation efficiency. The steep slope 5 means that most protoglobular clouds achieve 

' too low a star formation efficiency to give rise to bound star clusters following gas 

^ \ removal. As a result, most newly formed stars are scattered into the field soon after 

their formation. Gas removal during star formation in massive clouds is thus likely the 
prime cause of the predominance of field stars in the Galactic halo. The shape of the 
present-day cluster mass function depends weakly on the underlying distribution of 
the star formation efficiency. Finally, we show that a Gaussian mass function for the 
protoglobular clouds with a mean logmg ~ 6.1 — 6.2 and a standard deviation a < 0.4 
provides results very similar to those resulting from a truncated power-law cloud mass 
spectrum, that is, the distribution function of masses of protoglobular clouds influences 
only weakly the shape of the resulting globular star cluster initial mass function. The 
gas removal process and the protoglobular cloud mass-scale dominate the relevant 
physics. 

Key words: globular clusters: general - Galaxy: halo - Galaxy: formation 



o 
o 



1 INTRODUCTION 

Globular clusters are dense spherical gravitationally-bound 
stellar clusters. By virtue of the age of the oldest clusters 
(~ 13 Gyr), they are invaluable probes into the earliest evo- 
lutionary stages of their host galaxy. In that context how- 
ever, the original properties of (systems of) globular clus- 



ters will have been modified by the effects of a Hubble-time 
of evolution in their galactic environment. It is essential to 
disentangle formation from evolution. How the present-day 
globular cluster mass distribution in a large galaxy compares 
with the initial one constitutes one such striking example. 
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The cluster mass function ^, namely, the cluster number 
per constant logarithmic cluster mass interval dN/dlogm, 
which is proportional to the number of objects per magni- 
tude unit, constitutes a primary characteristic of any globu- 
lar cluster system hosted by a massive galaxy. Intriguingly, it 
shows only a weak dependence on the size, the morpholog- 
ical type or the environment of the host galaxy (Ashman 
& Zepf 1998, Harris 1999). This universal globular clus- 
ter mass function is usually fitted with a Gaussian with 
a mean of logm ~ 5.2 — 5.3 and a standard deviation of 
o"iogm — 0.5 — 0.6. The underlying mass spectrum (i.e., the 
number of objects per linear mass interval) is well described 
by a two-index power-law, with exponents ~ —2 and ~ —0.2 
above and below ~ 2 x 10 ' Mq, respectively (McLaughlin 
1994). The peak of the Gaussian function in fact coincides 
with the cluster mass at which the slope of the mass spec- 
trum changes. 

Individual globular clusters and globular cluster sys- 
tems having evolved over a Hubble-time in their galactic 
environment, the globular cluster initial mass function has 
remained model-dependent, with two competing hypothe- 
ses. It may have been a featureless power-law with a slope 
of ~ —1, the Gaussian function characteristic of old globular 
cluster populations then resulting from evolutionary effects, 
predominantly the preferential removal of the more vulnera- 
ble low- mass clusters (Fall & Zhang 2001). In that case, the 
cluster mass at the turnover of the present-day mass function 
depends on the age of the cluster system as well as on the 
cluster disruption time-scale in the galaxy of relevance. The 
older the cluster system and/or the shorter the disruption 
time-scale, the higher the cluster mass at the mass function 
turnover. In this situation, Vesperini (2001) shows that it is 
not straightforward to produce (almost) universal globular 
cluster mciss functions in very different types of galEixies if 
starting from an initial featureless power-law. Alternatively, 
Vesperini (1998) demonstrates that the present-day cluster 
mass function represents an equilibrium state "able to pre- 
serve its initial shape and parameters for one Hubble-time 
through a subtle balance between disruption of clusters and 
evolution of the masses of the surviving ones" , even though a 
significant fraction of the initial cluster population has been 
destroyed. That is, the initial cluster mass function may also 
be a Gaussian similar to that today. Should that be the 
case, the observed Gaussian shape of the cluster mass func- 
tion and its universality among galaxies are the preserved 
imprints of the globular cluster formation process. 

Parmentier & Gilmore (2005) and Vesperini et 
al. (2003) provide evidence for a Gaussian cluster initial 
mass function in the Galactic halo and in the giant ellip- 
tical M87, respectively. Theoretical support for a general 
bell-shaped cluster initial mass function has been missing so 
far however. 

The cluster initial mass function has often been as- 
sumed to mirror the mass function of the cluster gaseous 
progenitors (McLaughlin & Pudritz 1996, Elmegreen & Fal- 



garone 1996, Elmegreen & Efremov 1997). However, this is 
true only if all the cluster forming clouds turn the gas they 
are made of into bound star clusters with roughly the same 
efficiency. At its earliest stages, a protostellar cluster is em- 
bedded in the left-over star forming gas. That residual gas 
is blown away once the newly formed massive stars explode 
as supernovae. As a result of the corresponding weakening 
of its gravitational potential, the protocluster cam emerge 
unbound from the violent relaxation phase which follows. 
By means of analytical computations, Hills (1980) predicted 
that, in the case of instantaneous gas removal ^, initially 
virialized systems eventually dissolve for star forming effi- 
ciencies (i.e., the mass fraction of gas ending up in stars, 
hereafter e) smaller than 50%. That is, star forming clouds 
must be better than 50 % efficient in converting gas into 
stars in order to produce bound stellar clusters. The limited 
variations in e (i.e., less than a factor of 2) may then guar- 
antee that the initial mass function of the clusters is that of 
their parent clouds, shifted downwards by a factor of < 2. 

Using Af-body simulations however, Lada, Margulis & 
Deaxborn (1984) revisited that issue and pointed out that a 
system which becomes globally unbound due to a less than 
50 % efficiency is not necessarily completely disrupted and 
may retain a core of bound stars. Conversely, even a star 
forming efficiency larger than 50 % does not prevent the pro- 
tocluster from losing a fraction of its stars. This is due to 
the fact that the protocluster stars are characterized by a 
velocity distribution. The stars in the low- velocity tail of the 
distribution tend to survive as a gravitationally bound en- 
tity even if e < 0.5, while high-velocity stars escape even if 
e > 0.5. As a result, the initial mass of a stellar cluster is 
not determined by the star formation efficiency only. It de- 
pends on the mass fraction of the cluster parent cloud which 
is turned into stars which remain bound after the dispersal 
of the gaseous component. 

Modelling the dynamical evolution of gas-embedded 
clusters, that is, how they respond to supernova activ- 
ity, has since then attracted considerable attention (Ver- 
schueren 1989, Goodwin 1997, Geyer & Burkert 2001, Boily 
& Kroupa 2003, Fellhauer & Kroupa 2005). Surprisingly 
however, whether the early evolution of gas-embedded pro- 
toclusters influences the cluster initial mass function, that is, 
the mass function of the emerging gas-free bound groups of 
stars, has remained poorly explored (Kroupa & Boily 2002). 
Since the initial mass of a star cluster depends on gas re- 
moval, gas mass loss is actually worth being put to the test 
as a possible controlling mechanism of the globular cluster 
initial mass function. 

The outline of the paper is as follows. In section 2, we 
present the results of Monte-Carlo simulations highlighting 
the impact of gas removal upon the shape of the globular 
cluster initial mass function with respect to that of their 
gaseous progenitors. We also discuss how variations in the 
different model parameters affect the cluster initial mass 
function. In Sect. 3, we apply our model to the specific case 



^ In what follows, we adopt the nomenclature of McLaughlin & 
Pudritz (1996). We call mass spectrum the number of objects per 
linear mass Interval dN/dm, while we refer to the mass function 
to describe the number of objects per logarithmic mass interval 
dAT/dlog m. 



The gas removal takes place on a time-scale Tgr shorter than 
a protocluster crossing-time Tcross, implying that the stars do 
not have time to adjust to the new gravitational potential. Their 
velocity dispersions before and immediately after gas loss are thus 
considered to be the same. 
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of the Galactic stellar halo. Finally, we present our conclu- 
sions in Sec. 4. 



2 IMPACT OF GAS REMOVAL ON THE 
CLUSTER INITIAL MASS FUNCTION 

The formation of a stellar cluster is terminated when the 

newly formed massive stars go supernovac and blow away 
the gas left-over by the star formation process. Following 
the dispersal of that residual gas, the stars suddenly find 
themselves in a shallower gravitational potential, entailing 
either the escape of some of them, or even the complete 
destruction of the protocluster. If the gas is removed ex- 
plosively (i.e. Tgr « Tcross), the star formation efficiency 
must be larger than a threshold value eth = 33 % (see below 
and Fig. 1) for the protocluster to retain a bound core of 
stars. The mass fraction Fhound of stars remaining bound 
after gas removal ranges from zero (when the efficiency e 
is at its threshold value, or lower, i.e. e < 0.33) up to unity 
(when e < 1, so that gas removal is just a small perturbation 
of the stellar system). Considering the case of initially viri- 
alized gas-embedded protoclusters, various studies (Lada et 
al. 1984, Goycr & Burkort 2001, Boily & Kroupa 2003, FcU- 
hauer & Kroupa 2005) have led to fairly consistent results 
regarding the Fbound vs. e relation (see the plain symbols in 
Fig. 1). The knowledge of this relation enables us to relate 
the initial mass rriinit of a gas-free bound star cluster to the 
mass rricioud of its gaseous progenitor, namely: 



'^init — FliQund X € X TTldoud 



(1) 



As a result of the large variations in the bound star for- 
mation efficiency Fio^nd x e, an assumed simple mapping 
between the mass function of the cluster gaseous precursors 
on the one hand and the initial mass function of the clusters 
on the other hand can no longer be taken for granted. We 
now investigate how the cluster initial mass function differs 
with respect to the protoglobulax cloud mass function as a 
result of gas removal. 



2.1 From a truncated power- law mass function to 
a Gaussian mciss function 

In what follows, we assume that the protoglobular cloud 

mass spectrum obeys a power-law 



dN oc rn"drn , 



(2) 



with a varying between —2.5 and —1.5, as is observed for 
giant molecular clouds and their star forming cores in the 
Local Group of galaxies (e.g. Rosolowski 2005, see also sec- 
tion 3.1). A power-law mass spectrum may result from the 
coalescence of initially small equal-mass cores into a sys- 
tem of more massive objects with a wide range of masses 
(McLaughlin & Pudritz 1996). Alternatively, Elmegreen & 
Falgarone (1996) and Elmegreen & Efremov (1997) suggest 
that a power- law mass spectrum with a < —2 for the clus- 
ter gaseous progenitors is an imprint of the fractal structure 
of the star forming gas. As for the lower and upper limits 
of the cloud mass range, we adopt miow = 4 x W^Mq and 
rriup = IO'^Mq, that is, the Jeans mass range (see section 




Figure 1 . Relations between the fraction F(,o„„ci of stars remain- 
ing bound to the protocluster after gas removal and the star 
formation efRciency e achieved by the gaseous progenitor. The 
plain/open symbols correspond to the case of rapid/slow gas re- 
moval (i.c, Tgr << Tcross Or Tgr >> Tcross)- Data are from Lada 
et al. (1984), Gcycr & Burkcrt (2001), Boily & Kroupa (2003) 
and Fellhauer & Kroupa (2005) (respectively quoted as LMD84, 
GBOl, BK03 and FK05). The soUd and dotted lines depict the 
relations used in our simulations. 



3.5). In the next section, we will explore how the mass func- 
tion of the newly-born star clusters depends on these mass 
limits. 

Star forming regions arc characterized by a range in 
their respective star formation efficiency e, so that the pro- 
toglobular cloud mass spectrum is convolved with an e prob- 
ability distribution function, which we describe by a decreas- 
ing power-law of slope 6 and core rc, that is: 

ci I i H \ +C4 . (3) 



, , dN 
Pie) = 



1 + ^ 

Vc 



The two parameters ci and C4 arc determined so as to sat- 
isfy the two following constraints; (1) the integration of the 
probability distribution over the range e =[0,1] is unity, and 
the probability p(e) is zero when e = 1. 

The formation of a bound star cluster requires its 
gaseous progenitor to achieve e > eth (where "th" stands 
for "threshold"), i.e., the local star formation efficiency must 
be greater than ~ 0.3 — 0.4. On the scale of a galaxy, star 
formation proceeds inefficiently, so the global star formation 
efficiency may be of order a few per cent only. The core 
Vc and the slope S of the efficiency distribution p(e) arc 
thus bounded so that the mean star formation efficiency e, 
namely, the mass fraction of gas converted into stars for an 
entire system of protoglobular clouds, is one per cent: 



Jo 



£ diV(e) = 0.01 



(4) 



For a given system of protoglobular clouds, the slope S 
determines the fraction f>th of clouds achieving a star for- 
mation efficiency larger than the threshold eth and, there- 
fore, the initial size of the cluster system. Combined with 
the Fbound vs. e relation, it also determines the mean bound 
star formation efficiency e^. 



f 



(e X Fbound)dN{e) , 



(5) 
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Table 1. Fraction fyth of piotoglobular clouds giving rise to 
bound stellar clusters and mean bound star formation efficiency 
C5 for various distributions p(e) and thresholds c^h of the star for- 
mation efficiency, that is, for various triples (6, Tc, tth)- If <5 7^ 0, 
the core Tc and the slope 5 are bounded so that the mean star for- 
mation efficiency e across the whole galaxy is 0.01. The threshold 
etfi depends on the gas removal timescale Tgr (see section 2.2.3) 
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that is, the total mass fraction of gas converted into stars 
residing in bound systems after gas removal. To illustrate 
this with specific examples, Table 1 lists the values of f>th 
and for different triples {6, rc, eth) (see section 2.2.3 for 
the case of non-explosive gas removal leading to eth signifi- 
cantly less than 0.33). 

In our model, each protoglobular cloud is randomly as- 
signed a mass mdoud and a star formation efficiency e, both 
parameters being drawn independently from their respec- 
tive distribution, i.e. we assume no correlation between the 
mass of a cluster gaseous progenitor and its star formation 
efficiency. 

What the star formation efficiency of a star forming re- 
gion depends on is still ill-known. Elmegreen & Efremov 
(1997) suggest that, for any given ambient pressure, there 
is a trend of increasing star formation efficiency for more 
massive clouds (bottom panel of their Fig. 4). Consider- 
ing a large galaxy, such a relation might be useful in any 
region which is spatially limited enough so that the ambi- 
ent pressure does not vary significantly, that is, all clouds 
of that region are bounded by roughly the same pressure. 
However, our study is concerned with the formation of an 
entire star cluster system, spanning several tens of kpc in 
size and, therefore, spatial variations in the ambient pres- 
sure of the protoglobular clouds should be accounted for as 
well. Actually, Elmegreen & Efremov 's (1997) model also 
predicts that, for any given cloud mass, a higher ambient 
pressure promotes a larger star formation efficiency. As a 
result, two equal-mass protoglobular clouds may undergo 
markedly different star formation efficiencies due to differ- 
ent spatial locations in the protogalaxy, i.e., different exter- 
nal pressures. Clearly, our model cannot rely on a simple 
one-to-one radoud-e relation. 

What is more, the very existence of a monotonic mdoud- 
e relation, even within a spatially limited environment, is 
rather uncertain. Overplotting the star formation efficien- 
cies inferred for a few star forming regions of the solar neigh- 
bourhood (Lada & Lada 2003, their Table 2) on the mdoud-e 
relations predicted by Elmegreen & Efremov (1997), we note 
a large scatter and the inability of the mdoud-e relation cor- 
responding to the solar neighbourhood pressure to account 
for all data. Moreover, the presumed steady rise of the star 
formation efficiency with the cloud mass (for a given pres- 



sure) stems from an increase of the cloud specific binding 
energy, hence from an increase of the cloud velocity disper- 
sion (section 4 of Elmegreen & Efremov 1997). How the star 
formation efficiency responds to the latter is not accounted 
for, however. Higher gas velocity dispersions are promoted 
by larger turbulence and/or magnetic field pressures within 
the star forming cloud (e.g. Harris & Pudritz 1994), two 
processes known for hampering star formation. Specifically, 
Schmeja, Klessen & Proebrich 's (2005, their section 4.2) 
sinmlations demonstrate the existence of an inverse correla- 
tion of the star formation efficiency with the Mach number, 
that is, a higher velocity dispersion gas is less efficient at 
forming stars. Therefore, two counteracting effects arc at 
work: while the increase of the cloud specific binding energy 
raises the star formation efficiency, this gets lowered by the 
corresponding increase of the gas velocity dispersion, equiv- 
alently the growth of the magnetic fields and/or turbulence 
pervading the star forming gas. Even for a given ambient 
pressure, the star formation efficiency does not depend on 
the cloud mass monotonically. 

In the presence of significant variations in the ambient 
pressure, turbulence and magnetic field pressures from one 
star forming cloud to another, as may be expected when con- 
sidering a whole protogalaxy, the nidoud-e relation is likely 
reduced to a mere scatter, hence justifying our choice of a 
probability distribution function p{e) independent of any ex- 
plicit protoglobular cloud mass dependence. That is, given 
the lack of a clear physical prediction, and the observational 
uncertainties which we review in section 3.5 below, in this 
analysis we adopt a probabilistic approach, to investigate 
the extent to which it can be a viable option. 

Following the onset of supernova activity, the gas- 
embedded cluster gets exposed as its residual g£is is removed. 
Not only does the protocluster lose its gaseous component, 
it also loses a fraction of its initial stellar mass. We account 
for this phase by matching each efficiency value e to the cor- 
responding fraction Fbound of bound stars. The F^ound vs. e 
relation we are using is shown as the solid line in Fig. 1. It 
matches the results of Geyer & Burkert (2001) and Fellhauer 
& Kroupa (2005) and is characterized by eth ~ 0.33. 

Finally, the initial mass mi„it of globular clusters is de- 
rived following equation 1. The corresponding cluster initial 
mass functions, along with their protoglobular cloud mass 
functions, are shown in the top and middle panels of Fig. 2, 
assuming 5 = —4 and S = —2. It is worth emphasising that, 
although the simulations were started with truncated power- 
law cloud mass functions, the newly formed gas-free bound 
star clusters show bell-shaped mass functions. That is, lit- 
tle memory of the cloud mass function is retained as gas 
removal generates the low-mass regime of the cluster mass 
distribution. Therefore, gas removal is supported as a prime 
candidate mechanism responsible for generating bell-shaped 
globular cluster initial mass functions and, as such, may be 
central to our understanding of the origin of the bell-shaped 
universal globular cluster mass function. 

2.2 A possible origin for the universal location of 
the globuleir cluster mass function turnover 

In this section, we investigate how the derived cluster initial 
mass functions respond to model parameter variations. The 
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results are illustrated in the panels of Figs. 2-6. Unless oth- 
erwise stated on the panels, we adopt miow = 4 x 10 ' M©, 
ruup = 10^ M©, 6 = -2, eth = 0.33 and e = 0.01. For each 
simulation, Table 2 lists the values of the input parameters 
along with the characteristics of the generated cluster initial 
mass functions. Specifically, we present the logarithm of the 
cluster mass at the turnover logmTO, the mean logaritmic 
cluster mass logm, the standard deviation a, the skewness 
and the curtosis ^ of the logrri distribution. The cluster ini- 
tial mass functions arc not strictly Gaussian, as shown by the 
top and middle panels of Fig. 2-3, where we have overlaid 
the high mass regime of each mass function with a Gaus- 
sian. Our model predicts a mimbcr of low-mass (i.e. with 
a mass smaller than that at the turnover) clusters which is 
larger than if the cluster mass were drawn from a Gaussian 
distribution. As a result, the skewness of the cluster mass 
distributions is negative. Regarding the width of the mass 
function, Table 2 shows that the standard deviation a re- 
mains roughly constant regardless of the input parameter 
values. It is on the order of 0.6 and, therefore, consistent 
with the width of observed globular cluster mass functions. 

As emphasized in the introduction, if the globular clus- 
ter initial mass function is actually a bell-shape similar to 
that observed today, then the origin of the almost universal 
cluster mass at the turnover is locked into the cluster for- 
mation process. In the frame of our model, we now explore 
what the turnover location depends on. 

2.2.1 The spectral index a of the cloud mass spectrum 

We firstly address the eff'cct of varying the spectral index 
Q of the protoglobular cloud mass spectrum. We consider 
three values of a: —2.5, —2, —1.5. As a result of the greater 
fraction of high-mass clouds in the case of a shallower cloud 
mass spectrum, the cluster mass at the turnover may be 
larger if a = —1.5 than if a = —2.5. Such an effect actually 
shows up in case of a steep star formation efficiency distribu- 
tion {5 = —4), although the effect remains moderate, with 
[logmTo]a=-i.5 — [logmTo]Q!=-2.5 < 0.2. For shallower star 
formation efficiency distributions {S = —2 or 5 = 0), the 
effect is negligible (see Fig. 2). Memory of a steeper cloud 
mass spectrum is sometimes retained as a more pronounced 
skewness of the cluster initial mass function. 

2.2.2 The slope 5 of the star formation efficiency 

distribution p(e) 

Star forming regions do show variations in their respective 
star formation efficiency (Lada & Lada 2003) . In our model, 
this is accounted for by the functional form p(e), which de- 
scribes the probability distribution of the star formation ef- 
ficiency e. This is likely the most ill-dctcrmincd ingredient 
of our model. Actually, very few e mcEisurements exist for 
star forming regions in the Galactic disc since these require 
estimates of both the gaseous mass and the stellar mass. 
We know however that on a galactic scale, equivalent here 

^ We remind the reader that the skewness of a distribution char- 
acterizes its degree of asymetry while the curtosis measures its 
relative peakedness or flatness with respect to the null value of a 
Gaussian (Press et al. 1992). 



to the scale of a system of protoglobular clouds, star for- 
mation proceeds with a global efficiency of a few per cent 
only. Additionally, the vast majority of star forming regions 
give rise to unbound stellar groups, a process sometimes 
referred to as "infant mortality". In the Galactic disc for in- 
stance, the birthrate of embedded clusters in molecular cloud 
cores is observed to be extremely high compared to the birth 
rate of classical open clusters, thus suggesting that only a 
small fraction (a few per cent at most) of embedded clusters 
emerge from their natal cores as bound clusters. As quoted 
by Lada & Lada (2003), this high infant mortality results 
from the low to modest star formation efficiency and rapid 
gas dispersal that characterize their birth, that is, e < eth 
for most star forming regions. Observations of violent star 
forming regions in interacting and merging galaxies lead to 
similar conclusions. For instance. Fall, Chandar & Whitmore 
(2005) show that the age distribution of star clusters in the 
Antennae galaxies (NGC4038/39) declines approximately as 
dN/dr oc r^^ over the range 10® < r < 10'^ yr. They inter- 
pret this steep decline as evidence for a high rate of infant 
mortality, that is, most of the young clusters are not tightly 
gravitationally bound and are disrupted shortly after they 
form by the energy and momentum input from young mas- 
sive stars to the residual star forming gas, and massive star 
mass loss. 



Therefore, even though the overall distribution p(e) re- 
mains poorly determined, it actually makes sense to describe 
it as a decreasing power-law of the efficiency e with a slope 
steep enough so that only a small fraction of the protoglob- 
ular clouds convert gas into stars with an efficiency e larger 
than the threshold tth- We consider two values for the slope 
of p(e): 5 = —4 and S = —2. At the same time, we retain the 
constraint of a global star formation efficiency e of one per 
cent and wc determine rc accordingly (see Table 1). Addi- 
tionally, we investigate what happens if the star formation 
efficiency is uniformly distributed (i.e. p{e) is flat). In that 
case also, bell-shaped cluster initial mass functions emerge, 
although with a much sharper peak. The corresponding clus- 
ter initial mass functions arc shown in the three panels of 
Fig. 2 with, from top to bottom, 5 = —4, S = —2 and 6 = 0. 

As expected, the steeper the distribution p{e), the 
smaller the cluster mass at the turnover. The dependence 
of the turnover location on 5 remains limited however, with 
the steep {5 = —4) and flat {6 = 0) p(e) distributions 
leading to a difference of at most 0.4 in logm. The clus- 
ter mass at the turnover is even more robust in case of 
a shallow cloud mass spectrum (i.e. if —2 < a < —1.5, 
[logmTo]5=o — [logmTo]5=-4 ^ 0.3). If the mass spectrum 
of the protoglobular clouds were similar to that of the cluster 
forming cores of molecular clouds in the present-day Galac- 
tic disc, for which a ~ —1.7 (Lada & Lada 2003), then the 
uncertainties in the predicted turnover location arising from 
our misknowledge of p(e) are not larger than the uncertain- 
ties in the globular cluster mass-to-light ratio, which is on 
the order of a factor of 2. Varjdug the slope 5 of p(e) mostly 
affects the curtosis of the cluster initial mass function in 
the sense that a shallower efficiency distribution increases 
its peakedness. 
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Figure 2. Truncated power-law protoglobulaj: cloud mass func- 
tions and resulting bell-shaped cluster initial mass functions in 

the case of instantaneous gas removal (tjj, = 0.33, solid line in 
Fig. 1). The cloud mass range is 4 X 10^ Mq < rridoud < lO'^ Mq. 
The solid, dotted and daslied-dottcd lines correspond to a = 
—2.5, —2, —1.5, respectively, where a is the spectral index of the 
cloud mass spectrum. The bottom, middle and top panels cor- 
respond to three distinct slopes S of the probability distribution 
p{e) for the star formation efficiency, respectively, 5 = —4, —2 and 
0. In the top and middle panels, Gaussian curves overlaid on each 
cluster initial mass function show that our model predicts more 
low-mass clusters than does a pure Gaussian mass function. 



2.2.3 The star formation efficiency threshold eth 

The ultimate fate of an embedded cluster (i.e. whether it will 
survive as a bound stellar cluster or not) depends on whether 
the star formation cfBcicncy e achieved by the gaseous pro- 
cursor is larger than the threshold ttn or not. The efficiency 
threshold itself depends on the timescale Tgr for remov- 
ing the gas out of the protocluster. If the gas is removed 
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Figure 3. Same as Fig. 2, but in case of slow gas removal (e^/j 
0.15, dotted line in Fig. 1). 



instantaneously (i.e., Tgr « Tcross), we get eth — 0.33 
(see the plain symbols in Fig. 1). For slower gas removal 
[Tgr > Tcross) however, the efficiency threshold gets smaller 
(sec the open symbols in Fig. 1, data taken from Geyer & 
Burkcrt 2001) because the stars can now adjust adiabati- 
cally to the new gravitational potential they sit iu and ex- 
pand to a new state of virial equilibrium, even for e < 0.33. 

If Tgr ~ 10 X Tcross, eth — 0.15. 

Although protoclustcrs containing O stars arc expected 
to remove any residual star forming gas on a timescale 
shorter than Tcross (Goyer & Burkert 2001, Lada & Lada 
2003, Kroupa 2005), it is nevertheless interesting to inves- 
tigate whether a slow gas dispersal affects our results sig- 
nificantly. This situation mimics a weakly bound cluster 
where mass loss from the most massive stars on a stel- 



lar evolution timescale (~ 10 yr) can be significant. Us- 
ing the dashed curve in Fig. 1 as the Fbound vs. e rela^ 
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Figure 4. Same as the middle panel of Fig. 2, but for three 
different lower limits miaui oi the cloud mass range. The mass 
at the turnover of the cluster initial mass function sensitively 
depends on mjou,. 



tion (i.e. considering Tgr — IOtctoss and tth = 0.15 instead 
of Tgr « Tcross and eth = 0.33), we obtain the cluster 
initial mass functions of Fig. 3, each panel differing from 
its counterpart in Fig. 2 by the star formation efficiency 
threshold only. With respect to the case of instantaneous 
gas expulsion, any vahie of Ftound is now coupled with a 
lower value of e and we therefore expect a downward shift of 
the turnover location. This remains moderate however, with 

[logmTo]etft=0.33 - [logmTo]et;.=0.16 < 0.25. 



2.2.4 The lower limit niiow of the cloud mass range 

Unlike a and S, the lower limit of the protoglobular cloud 
mass range is a prime controlling parameter of the turnover 
location, as highlighted in Fig. 4. For instance, a four times 



smaller lower mass limit (i.e. 10^ Mq instead of 4 x lO'' Mq) 
results in a turnover shifted by —0.6 in logm (compare top 
and middle panels of Fig. 4). The turnover of the cluster ini- 
tial mass function thus tracks the lower mass limit of the pro- 
genitor clouds, with the cluster mass rriro at the turnover 
being of order the lower limit miow of the protoglobular 
cloud mass spectrum. Specifically, the difference miow—iriTO 
depends on the slope a of the cloud mass spectrum, on the 
slope 5 of the star formation efBciency distribution p(e) and 
on the efficiency threshold eth. A steeper cloud mass spec- 
trum, a steeper star formation efficiency distribution and/or 
a longer gas removal time-scale increase the offset between 
the cluster mass at the turnover and the lower cloud mass 
limit. 



2.2.5 The upper limit of the cloud mass range 

We now explore the effect of varying the uppper limit m.up of 
the protoglobular cloud mass range. In Fig. 5, three values of 
rriup are considered, from top to bottom: 10^ Mq , 3 x 10^ Mq 
and 10® Mq. The turnover location is practically unaffected, 
the shift in logm being < 0.05 when m^p is decreased by a 
factor of ten. The upper cloud mass limit however influ- 
ences markedly the high-mass regime of the cluster initial 
mass function, an effect which we will investigate in detail 
in section 3.2. 



2.2.6 The averaged star formation efficiency e 

Giant molecular clouds in the present-day Galactic disc turn 
their gas into stars with a global efficiency in the range be- 
tween 1 to 5 per cent (Lada &: Lada 2003). High star for- 
mation efficiencies are achieved only locally, in some dense 
cores. All the simulations performed up to now assume a 
global efficiency e = 0.01. We now check whether a five times 
larger value e = 0.05 affects the shape of the cluster initial 
mass function significantly. Figure 6 is the counterpart of the 
middle panel of Fig. 2. Comparison of both panels demon- 
strates the robustness of the turnover location. In fact, the 
increase in e is accounted for by a core radius of the efficiency 
distribution p{e) ten times greater (i.e, rc — 0.02) than that 
derived if e = 0.01. The efficiency distribution p(e) thus re- 
mains almost unaffected above the threshold eth and so does 
the cluster initial mass function. 

It follows from the above simulations that a bell-shaped 
cluster initial mass function arises from gas removal, pro- 
vided that the power-law protoglobular cloud mass spectrum 
shows a lower-mass limit. Assuming that the star formation 
efficiency is independent of the protoglobular cloud mass 
(see section 2.1), the turnover location depends weakly only 
on all but one of the model parameters, namely, the lower 
mass-limit of the protoglobular cloud mass spectrum. To 
first order, the cloud lower mass limit dictates the turnover 
location, that is, the lower miom, the smaller the cluster 
mass rriTO at the turnover. As an extreme case, a lower 
cloud mass limit of, say, miow = 100 Mq leads to a turnover 
located at mro ^ 100 Mq. If the cluster detection limit is 
higher than mro, then the observed cluster mass spectrum 
is a power-law whose slope is similar to that of the gaseous 
progenitors, by virtue of the mass- independent p(e) distribu- 
tion. In section 3.5, we will discuss further the importance of 
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Figure 5. Same as the middle panel of Fig. 2, but for three 
different upper limits rriup of the cloud mass range. The turnover 
location is unaffected by m^p variations. 



an appropriate protoglobular cloud mass-scale with respect 
to explaining the universal Gaussian globular cluster mass 
function. 



3 AN APPLICATION TO THE GALACTIC 
STELLAR HALO 

3.1 Fitting the halo cluster mass function 

In this section, we explore whether our model can success- 
fully reproduce the observed mass function of the Galactic 
halo globular cluster system. In order to evolve up to an age 
of 13 Gyr the cluster initial mass functions derived in the 
previous section, we make use of Baumgardt & Makino 's 
(2003) equation 12, which they derived by fitting the results 
of a large set of iV-body simulations taking into account 
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Figure 6. Same as the middle panel of Fig. 2, but with the core 
rc of the star formation efficiency distribution p(e) adjusted so 
that the global mean efficiency e is 0.05 instead of 0.01. 



the effects of stellar evolution, of two-body relaxation and 
of cluster tidal truncation. This analytical formula supplies 
at any time t the mass of a star cluster with initial mass 
rrii which is moving along a cirular orbit at a galactocentric 
distance D: 



m{t) 



= 0.70 1 



t 

idis 



(6) 



This equation assumes that mass loss due to massive-star 
stellar evolution takes place immediately after cluster forma- 
tion and amounts to 30 per cent of the cluster initial mass. 
Later on, globular cluster mass decreases linearly with time 
due to the combination of internal (two-body relaxation) 
and external (galactic tidal field) effects. The disruption 
time tdis of a cluster (equation 10 of Baumgardt & Makino 
2003) depends on the circular velocity Vc of the host galaxy, 
on the cluster initial mass rrii and on the cluster galactocen- 
tric distance D. In what follows, we adopt Vc = 220km.s~^ 
and we distribute the globular clusters in space following 
an initial number density profile scaling eis D~^'^ (Parmen- 
tier & Gilmore 2005). Additionally, we account for the non- 
circularity of the globular cluster orbits by halving the clus- 
ter disruption time-scale, which is equivalent to assuming 
a mean orbital eccentricity e = 0.5 (Baumgardt & Makino 
2003, their equation 10). Since the globular cluster mass 
function is an equilibrium shape (see section 3.3), that is, 
its shape is practically independent of the age t of the clus- 
ter system as well as of the exact cluster disruption time- 
scale tdis, the choice of that correction factor for tdis does 
not significantly affect the goodness of fit of the present-day 
globular cluster mass function. Only the fraction Fn of sur- 
viving clusters and the ratio Fm of the final to the initial 
mass in globular clusters are affected (see Table 3) . In order 
to take into account dynamical friction, clusters whose time- 
scale of orbital decay is smaller than t are removed from the 
cluster system at that time (Binney & Tremaine 1987). 

In this study, we focus on the mass distribution of 
the first generation globular clusters which formed within 
the gravitational potential well of the Galaxy. Hence, we 
do not consider the more metal-rich, presumably second- 
generation, bulge/disc globular clusters (Zinn 1985). The 
halo cluster system itself hosts two distinct populations of 
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Table 2. Characteristics of the initial mass functions of gas-free bound star clusters as derived utilising the model of equation 1 for 
various sets of model parameters. The first column gives the star formation efficiency threshold ct^ required for a protocluster to retain 
a bound core of stars (equivalently, the gas removal time-scale, see section 2.2.3). The slope 5 and the core Tc of the star formation 
efBciency distribution p{e) are listed in the second and third columns. Both parameters are adjusted so that the mean stax formation 
efficiency e is 0.01 (except in the case of a flat distribution p{e) for which 5 = 0). The next three columns give the lower limit mi„^, the 
upper limit ruup and the slope a of the protoglobular cloud mass spectrum. The last five columns list the logarithm of the cluster mass 
at the turnover logrnxo of the cluster initial mass function, its average value logm, standard deviation <t, skewness and curtosis. 
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clusters, the so-called Old Halo and Younger Halo (Van den 
Bergh 1993, Zinn 1993). The Old Halo globular clusters 
might have been formed 'in situ'. In contrast. Younger Halo 
globular clusters arc suspected of having been accreted. Re- 
gardless of their formation history, the Old Halo globular 
clusters form a coherent and well-defined group, well-suited 
to an analysis of their properties (see Table 1 in Parmen- 



ticr & Grebcl 2005 for the definition of our sample). Our 
source for the absolute visual magnitudes M„ of the clusters 
is the McMaster database compiled and maintained by Har- 
ris (1996, updated February 2003). Tlu;st; are converted into 
luminous mass estimates on the assumption of a constant 
cluster integrated mass-to-light ratio M/L^ of 2.35 (i.e., the 
average of the mass-to-light ratios of the halo clusters for 
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which Pryor & Meylan (1993) obtained dynamical mass es- 
timates) . 

Having detailed the cluster evolutionary model and the 
cluster sample against which we will compare the outcomes 
of our simulations, we now turn to the issue of deriving a 
cluster initial mass function which, when evolved to an age 
of 13 Gyr, reproduces the observed Old Halo cluster mass 
function. The shape of the cluster initial mass function de- 
pends on the following: 

(i) the truncated protoglobular cloud mass spectrum, de- 
fined by its slope a and its lower and upper mass limits miow 
and niup, 

(a) the star formation efficiency probability distribution 
p{c), defined by its slope 5 and its core Tc (equation 3), 
(in) the fraction of stars remaining bound to the protoclus- 
ter after gas removal, as defined by the Fbound vs. e diagram 
(Fig. 1). 

In what follows, we assume that any left-over star form- 
ing gas in each protocluster is instantaneously removed, and 
we therefore consider the solid line in Fig. 1 as the Fhound 
vs. e relation. For the slope of the cloud mass spectrum, 
we adopt a = —1.7, which is the spectral index reported 
by surveys of giant molecular clouds in a sample of galax- 
ies in the Local Group, namely, the outer disc of the Milky 
Way, the Large Magellanic Cloud, M31, and ICIO (Blitz et 
al. 2006, their Table 3). Although, shallower (a = -1.5) 
and steeper (a = —2.5) spectral indices are reported for the 
inner disc of the Milky Way and the spiral galaxy M33, re- 
spectively, these axe likely ascribed to observational biaises. 
As noted by Rosolowski (2005), in the inner Milky Way, 
line-of-sight blending will make several less massive clouds 
appear as a single more massive one, shifting the index to- 
ward shallower values. On the other hand, the steep slope 
of the cloud mass distribution in M33 can be attributed to 
fitting a truncated power-law distribution above the mass 
cutoff. In the present-day Galactic disc, the same spectral 
index a = —1.7 also describes the mass spectrum of the 
dense cores of giant molecular clouds, to which embedded 
clusters are often physically associated (Lada & Lada 2003 
and references therin). The lower (miow) and upper (niup) 
limits of the protoglobular cloud mass spectrum are left as 
free parameters. 

As already quoted in section 2.2.2, the star formation 
efficiency distribution p{e) is poorly determined. Moreover, 
even if it were reasonably known for the present-day Galactic 
disc, there is no guarantee that this would equally-well de- 
scribe star formation in the Galactic halo some 13 Gyr ago. 
Nevertheless, we attempt to constrain p{e) in the following 
way. We assume that the global star formation efficiency is 
on the order of a few per cent, thereby reflecting the inef- 
fectiveness of star formation of galactic scales. We note in 
passing that the ratio between the baryonic masses of the 
stellar halo and of the Galaxy is on that order of magni- 
tude. This may suggest that the stellar halo formed out of 
an amount of gas equivalent to the total mass in gas and 
stars in the present-day Milky Way. Thus, we adopt the fol- 
lowing mean star formation efficiency: 



In this equation, Mj^lj^ and M 



13 Gyr 
Halo 



are the initial and 



0.7" 



Halo 



present-day stellar masses of the Galactic halo, respectively, 
the factor 0.7~^ accounting for the corresponding stellar evo- 
lutionary mass loss. The stellar halo is processed out of an 
amount of gas M^"'' which, under the above mentioned as- 
sumption, is equivalent to the total mass M^^^^J^^^^ 
in stars and gas in the present-day Milky Way. We adopt 
a present-day halo mass M^^^^ of 10^ Mq (Freeman & 
Bland-Hawthorn 2002) and a present-day total mass of stars 
and gas in the Milky Way „. of 6 x lO" Mq, 

leading to a mean star formation efficiency e = 2.5 %. We 
acknowledge that this constitutes a very crude approxima- 
tion since equation 7 neglects the effects of later mergers and 
accretions onto the Milky Way. As we discuss below however 
(see Table 3 and the third panel of Fig. 8), the exact value 
of e does not significantly affect the goodness of fit to the 
Old Halo cluster mass function. 

As a second constraint on the star formation efficiency 
distribution p(e), we use the present-day mass fraction of 
Old Halo clusters in the Galactic halo, namely: 

rl3 Gyr _ ^V-itot 
/( 



IGC 



^13 Gyr 
^"Halo 



: 0.02 , 



(8) 



where is the present-day total mass in Old Halo globu- 

lar clusters, that is, 2 x lO'^ M© (Parmentier & Grebel 2005). 
The mass of the Galactic halo is actually dominated by field 
stars. In this model, these arise from (i) the disruption of 
protoclusters in the e < eth regime (infant mortality), (ii) 
the escape of stars out of protoclusters following gas removal 
when e > tth (infant weight- loss ) , and (in) the evaporation 
and disruption of globular clusters over a Hubble-time of 
dynamical evolution in the tidal field of the Milky Way. The 
present-day cluster mass fraction f^^^^ ^ the mass fraction 
Fm of surviving clusters at an age t = 13 Gyr, the mean star 
formation efficiency e and the mean bound star formation 
efficiency et are related through: 



f 13 Gyr 
JGC 



Fm X 66 
0.7 e 



(9) 



tb depends on the star formation efficiency distribution p(e) 
and on the Ftound vs. e relation, while e depends on the for- 
mer only (see equations 4 and 5). The top panel of Fig. 7 
shows the locii of points with identical mean star formation 
efficiencies e (thick curves) and the locii of points with iden- 
tical mean bound star formation efficiencies et (thin curves) 
in the 5 vs. rc diagram. This plot enables us to determine the 
function p(e), via its parameters d and rc, for any given cou- 
ple (e, eh). Yet, equation 9 shows that matching the distribu- 
tion p(e) and the present-day cluster mass fraction /, 



13 Gyr ■ 

GC 



jiflSGyr 

^"Gas+stars in MW 



(7) 



the halo still requires an estimate of Fm- The mass fraction 
Fm of surviving clusters depends on the initial mass function 
of the clusters, on their initial distribution in space, on their 
age and on the assumed cluster disruption time-scale tdis. 
The greater the contribution of the (more vulnerable) low- 
mass clusters, the more concentrated around the Galactic 
centre the cluster spatial distribution, the older the cluster 
system and/or the shorter the cluster disruption time-scale, 
the smaller the ratio Fm- We have detailed our assumptions 
regarding the last three quantities earlier in this section. We 
now need to determine a cluster initial mass function com- 
patible with equation 8. Since Fm is required to determine, 
in turn, (equation 9) , p(e) (via the top panel of Fig. 7) and 
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Figure 7. Top panel: Locii of points witli identical mean star 
formation efficiencies c (thick curves) and identical mean bound 
star formation efficiencies (thin curves) in the (rc,<5) plane, 
where S and rc are the slope and the core of the efficiency distri- 
bution p(e). Middle panel: Incomplete gamma functions Q of the 
fit of the evolved cluster mass function to the observed one (thick 
curves) and cluster mass fractions /gc ^-t ^S'^ of 13 Gyr as 
predicted by the model (thin curves) versus the lower and upper 
limits {niioiu, rfiup) of the cloud mass spectrum. The dotted area 
corresponds to simulations for which the mass at the turnover of 
the evolved cluster mass function is larger than 3 x 10^ Mq (see 
text for details). Bottom panel: Same as middle panel, but for 
a steeper p(e) function. This results in a greater rate of cluster 
infant mortality, thus lowering the present-day cluster mass frac- 
tion for any cloud mass range. The four-branch stars represent 
different cases presented in Fig. 8 and Table 3. The solid star is 
our adopted fiducial case (first line in Table 3). 



the cluster initial mass function itself (with the model de- 
scribed in section 2.1), we need to iterate. The middle panel 
of Fig. 7 summarizes the outcomes of a first set of simula- 
tions with {5, rc) — {-2A, 0.013). This corresponds to (e, it) = 
(0.025, 0.002). The x and y-axcs represent the lower and up- 
per limits of the protoglobular cloud mass spectrum, respec- 
tively. The thick curves are the locii of points with the same 
value for the incomplete gamma function Q. This quantifies 
the goodness of fit between the evolved cluster mass func- 
tion, as predicted by our model, and the Old Halo cluster 
mass function The dashed-dotted, solid and dotted lines 
correspond to Q=0.001, 0.1 and 0.3 respectively. The dot- 
ted region of the diagram, at log rniow ^ 5.9, corresponds to 
simulations for which the cluster mass at the turnover of the 
evolved mass function is larger than 3 x 10^ M©, that is, the 
turnover is located at too large a mass with respect to what 
is observed. This result seems to contradict the excellent 
goodness of fit (0.1 < Q < 0.3) in this region. The top panel 
of Fig. 8 highlights the origin of that discrepancy. The lower 
solid line is the modelled cluster mass function at an age of 
13 Gyr for our fiducial case (see below) and the plain squares 
depict the Old Halo cluster mass function. Our model over- 
estimates slightly the number of low-mass clusters. As il- 
lustrated in that panel, an increase in the lower limit miow 
of the cloud mass range shifts the cluster mass function to- 
ward larger mass. This improves the goodness of fit of the 
low-mass wing of the evolved cluster mass function and in- 
creases the incomplete gamma function Q accordingly, even 
though the cluster mass at the turnover gets too large (com- 
pare the lower solid line, for which miow = 6 x 10® M© , to the 
lower dashed-dotted line, corresponding to miow = 10^ M©). 
Coming back to the middle panel of Fig. 7, the two thin 
lines arc the locii of points with given present-day cluster 
mass fractions in the Galactic halo fQf?^^- The solid line 
(at the bottom left of the plot) and the dotted line corre- 
spond to fcfc^^^ = 0.02 and flfc?^^ = 0.03, respectively. 
The region of the plane {rniow,Tnup) where the evolved clus- 
ter mass function fits well the observed one (i.e. Q > 0.1, but 
outside the dotted area) coincides with a slightly too large 
cluster mass fraction, that is, f^(?^^ ^ 0.03. As shown in 
Table 1, a steeper slope 5 for the star formation efficiency 
distribution p(e) -while adjusting its core Tc so as to retain 
the same mean star formation efficiency e- results into a 
smaller mean bound star formation efficiency tb and, there- 
fore, a smaller cluster mass fraction f^(?^^ at an age of 
13 Gyr (sec equation 9). The bottom panel of Fig. 7 illus- 
trates the effect of steepening p(e), namely, 5 — —2.9 corre- 
sponding to €b = 0.001 (top panel of Fig. 7). As explained 
in section 2.2.2 (see also Figs. 2 and 3), the cluster mass 
function does not strongly depend on 5. Consequently, while 
the isofac curves arc moved upwards when <5 is decreased, 
the locus of points of any given value of Q is little affected. 
One can now see that our new choice of 5 enables us to fit 
the Old Halo cluster mass function (i.e. Q > 0.1) while, at 
the same time, accounting for the present-day cluster mass 



^ We remind the reader that: a Q value of 0.1 or larger indicates 
a satisfactory agreement between the model and the data; if Q > 
0.001, the fit may be acceptable if, e.g., the errors have been 
moderately underestimated; if Q < 0.001, the model can be called 
into question (Press et al. 1992). 
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fraction in the stellar halo ,fQ(?^^ — 0.02. In what follows, 
we adopt 5 = -2.9, = 0.025, miovj = 6 x lO'^M© and 
ruup = 5 X lO^M© as our fiducial case {Q = 0.1, depicted 
as the plain four-branch star in the bottom panel of Fig. 7) 
and wo explore the sensitivity of our results to variations in 
the input parameters {miow, mup, 5 and rc). For each sim- 
ulation. Table 3 lists the mean star formation efficiency e 
and the mean bound star formation efficiency et, the mass 
fraction f^(?^^ of clustered stars in the halo at an age of 
13 Gyr, the ratio Fm between the final total mass and the 
initial total mass in clusters, its number counterpart Fn 
(i.e. the fraction of surviving clusters) and the goodness of 
fit Q of the evolved cluster mass function to the observed 
one. The corresponding modelled cluster mass distributions, 
both initial and evolved to an age of 13 Gyr, are shown in 
the different panels of Figs. 8 and 9, together with the Old 
Halo cluster mass distribution. The solid lines correspond 
to the fiducial case. In all these panels, the initial cluster 
mass distributions have been corrected for the effect of stel- 
lar evolutionary mass loss, that is, all cluster masses have 
been decreased by 30% with respect to their actual initial 
value (see equation 6). This correction, indicated by a left- 
ward arrow of size 0.15 in logrn, enables us to estimate how 
the cluster mass function evolves with time due only to clus- 
ter dynamical evolution (i.e. evaporation and disruption). 

(i) Variations in miow 

As noted in section 2.2.4, the turnover location depends pri- 
marily on the lower limit of the cloud mass range. A smaller 
value of miow shifts the turnover downward, thereby raising 
the relative fraction of low-mass clusters (see top panel of 
Fig. 8, in which miow = 2 x 10^ 6 x 10^ and IO^Mq). The 
present-day mass fraction f^^^^ of clusters in the halo and 
the survival ratios Fm and Fn are therefore decreased with 
respect to the fiducial case. 

(ii) Variations in /q'^^^ 

Varying the slope 5 of the efficiency distribution p(e), while 
retaining the same mean star formation efficiency e, alters 
the fraction of protoglobular clouds giving rise to bound 
star clusters following gas removal, that is, those for which 
e > eth- As illustrated by the comparison between the mid- 
dle and bottom panels of Fig. 7, a smaller d value lessens 
the present-day mass fraction of clusters f^Q^^ through a 
lower mean bound star formation efficiency eE (see top panel 
and equation 9). Accordingly, the cluster initial mass func- 
tion reveals a larger fraction of low-mass clusters. Yet, at an 
age of 13 Gyr, the evolved cluster mass functions are undis- 
tinguishable (see the second panel of Fig. 8 for which 5 = 
-3.9,-2.9,-1.9, corresponding to f^c^'' ~ 0.01,0.020.05, 
respectively) . 

(iii) Variations in e 

We have already emphasized how poorly constrained is our 
estimate e = 2.5 % of the mean star formation efficiency. 
However, the third panel in Fig. 8 shows that a lower (1 %) 
or higher (5 %) value for e, while retaining the observational 
constraint f^c^^ = 2 %, hardly affects the shape of the 
evolved cluster mass function. 

(iv) Variations in tdis 

In the simulations described so far in the present section, 
we have assumed a mean eccentricity e = 0.5 for the cluster 
orbits. Equivalently, the disruption time-scale of any cluster 
is half (i.e. (1— e)) that of the same cluster on a circular orbit. 
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Figure 8. Initial/evolved (plain/open circles) cluster mass func- 
tions corresponding to the cases listed in Table 3. In each panel, 
the solid lines are the fiducial case. The arrow indicates that the 
initial mass functions have been shifted by -0.15, that is, corrected 
for stellaj: evolutionary mass loss. The top and bottom panels 
show the effect of varying the lower cloud mass limit ttiiow 

and 

the cluster disruption time-scale t^is, respectively. The middle 
panels display the results for different present-day cluster mass 
fractions f^Q^^ and different mean star formation efficiencies ? 
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Table 3. Impact of varying the protoglobular cloud mass range and the star formation efficiency distribution p(e) on tfie model outcomes. 
We successively explore how variations in the lower cloud mass limit miow, the upper cloud mass limit rriup, the slope S and core rc 
of p{e) affect the cluster mass fraction in the stellar halo 1^^^^ goodness of fit Q of the cluster mass function at an age of 

13 Gyr, as well as the ratio Fm between the final and initial masses in clusters and the fra<;tion Fj^ of surviving clusters, with respect 
to the fiducial case (first line of the Table). We also consider a cluster disruption time-scale t^j^ twice as high as for the fiducial case, 
corresponding to clusters on circular orbits. The mean (bound) star formation efficiency e (e^) is given for each couple (5, rc). 
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The bottom panel of Fig. 8 compares our fiducial case to the 
case for which all clusters are on circular orbits. It appears 
that the exact scahng of the cluster disruption time-scale 
(or, oquivalcntly, the exact ago of the cluster system) does 
not influence the shape of the cluster evolved mass function 
either. This is so because the observed globular cluster mass 
function is a state of dynamical equilibrium, its shape being 
preserved over the cluster system evolution (see section 3.3). 

Encounters between gas clouds and globular clusters dur- 
ing their early secular evolution, when they orbit through 
still gas-rich environments, constitute another source of 
uncertainty in the cluster disruption time-scale. Gieles et 
al. (2006) find that, when taking into account giant molecu- 
lar cloud encounters, the disruption time-scale of clusters in 
the solar neighbourhood is shortened by a factor ~ 3.5 com- 
pared to what is found in simulations taking into account the 
combined effect of Galactic tidal field and stellar evolution 
only. Again, this uncertainty is of little relevance to the evo- 
lution with time of the globular cluster mass function, since 
this represents an equilibrium state. This is exemplified by 
Fig. 3 in de Grijs, Parmentier & Lamers (2005) which shows 
that the evolved counterpart of a Gaussian cluster initial 
mass function similar to that inferred in this section is prac- 
tically the same, regardless of whether the cluster disruption 
time-scale is that predicted by Baumgardt & Makino (2003) 
or 10 times shorter. 

The main effect of a shortened cluster disruption time- 
scale is a lowering of the cluster survival rates Fm and Fn, 
equivalently, an increase of the initial number of clusters in 
order to explain any given present-day cluster population. 

(v) Variations in ruup 
How the upper limit niup of the cloud mass range affects the 
results is discussed in the next section. 

Examination of the panels of Figs. 8 and 9 reveals that, 
over the mass range m < 5 x 10* Mq , none of our mod- 
els fit the observed mass distribution in the sense that they 



all slightly overpredict the number of clusters. The discrep- 
ancy is far from severe, however, as it remains limited to 
less than a 2-a error. Whether this effect truly resides in 
our model or reflects the incompleteness of the cluster sam- 
ple is uncertain. A few globular clusters located on the other 
side of the Galactic bulge may have run undetected thus far 
(see, e.g.. Fig. 1 of Larsen 2006 which shows the asymmetri- 
cal distribution of observed globular clusters in the Galactic 
plane). Missing clusters tend to preferentially include low- 
mass ones, owing to the combination of their lower luminos- 
ity and lower concentration (see Fig. 4 of Larsen 2006 for the 
mass-concentration correlation). They may thus account for 
the discrepancy noticed in Figs. 8 and 9. This possibility 
gets strengthened by a comparison between the mass spec- 
trum of the Old Halo cluster system (plain squares in the top 
panel of Fig. 9) and that of the more populous M31 globular 
cluster system (see Fig. 1 of McLaughlin & Pudritz 1996). 
While the former tends to steadily decrease with decreasing 
cluster mass over the mass range m < 5 x 10* Mq, the latter 
remains flat, possibly even slightly increasing with decreas- 
ing cluster mass, thereby agreing with what is predicted by 
our model. 



3.2 Substructures in the high-mass regime of the 
cluster mass distribution 

Although the globular cluster mass at the turnover of the 

present-day Gaussian mass function is practically constant 
from one massive galaxy to another, the overall shape of 
the globular cluster mass function is not strictly univer- 
sal. For instance, McLaughlin & Pudritz (1996) note that 
the high-mass regime of the mass function of the globu- 
lar cluster systems of the Milky Way and M31 combined 
differs from that of the giant elliptical M87. These differ- 
ences are best revealed by the mass-weighted mass func- 
tion m^dN/dm oc mdiV/dlogm, which is derived by sum- 
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Figure 9. Initial/evolved (solid/open circles) cluster mass dis- 
tributions for three distinct upper cloud mass limits va^p (see 
Table 3). In each panel, the solid lines are the fiducial case. The 
arrow indicates that the initial mass functions have been shifted 
by -0.15, that is, corrected for stellar evolutionary mass loss. The 
Old Halo mass distribution is depicted as the plain squares. The 
top, middle and bottom panels show the mass spectra, the mass 
functions and the mass- weighted mass functions, respectively. 



ming individual cluster masses in evenly spaced logarith- 
mic mass bins. The mass-weighted mass function therefore 

describes how the total mass of the cluster system is dis- 
tributed between the high- and low-mass members of the 
cluster population. If the spectral index of the cluster mass 
spectrum diV/dm is a = —2, the mass-weighted mass func- 
tion is flat and any two identical-size logarithmic mass bins 
contribute equally to the total mass of the cluster system 
(e.g. the contributions to the system total mass of the clus- 
ters in the mass range 10^ < m < lO^M© and of those in 
the mass range 10^ < m < IO^Mq are the same). A shal- 
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Figure 10. Evolved cluster lognormal mass functions generated 
by our model (solid symbols) for various upper limits mup of 
the power-law cloud mass spectrum ( Toy panel, see Table 3) 
and various standard deviations a of the Gaussian cloud mass 
function (Bottom panel, see Table 4). The solid lines correspond 
to the functional form given by equation 10 and introduced by 
Burkert & Smith (2000) to describe globular cluster mass dis- 
tributions above the turnover of the mass function. Top panel: 
rUc = 0.5 X 1O®M0,O.9 x 1O6M0,2 x IO'^Mq- Bottom panel: 
rUc = l.l X 1O6M0,6 X IO^Mq. 



lower (steeper) mass spectrum implies that the total mass is 
mostly contained in the high-mass (low-mass) objects and 
that the mass- weighted mass function is increasing (decreas- 
ing) with cluster mass. 

McLaughlin & Pudritz (1996) note that there is a clus- 
ter mass iripeak, larger than that at the mass function 
turnover, at which the mass-weighted mass function peaks. 
That is, the spectral index of the mass spectrum steepens 
from Q > —2 to Of < —2 and, above the turnover, the glob- 
ular cluster mass spectrum cannot be accurately fitted with 
a single power-law. 

The mass-weighted mass function of the Old Halo, as 
well as those predicted by our model for various upper pro- 
toglobular cloud mass limits niup, are shown in the bottom 
panel of Fig. 9. The observed distribution reveals that the 
contribution of low- mass (< W^Mq) globular clusters to the 
total mass of the Old Halo cluster system is negligible. The 
absence of any peak in the mass-weighted mass function of 
the Old Halo likely arises from the limited size of the cluster 
sample (72 clusters). In contrast, the mass-weighted mass 
function of the globular cluster system of the giant ellipti- 
cal M87 shows a prominent peak (McLaughlin & Pudritz 
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1996, their Fig. 2). Owing to a cluster number about two 
orders of magnitude larger than that of the Old Halo, the 
M87 mass distribution is probed up to higher cluster mass 
and the substructures of the high-mass regime are therefore 
better highlighted. 

Interestingly enough, our model predicts a peak in the 
mass-weighted mass function. The bottom panel of Fig. 9 
shows that its location depends on the upper limit rUup of 
the protoglobular cloud mass range, the higher m„p, the 
higher nipeak- The dotted, solid and dashed-dottcd Hues cor- 
respond to rUup = 3 X 10®, 5 X 10®, 10 ""Mo, respectively. The 
corresponding mass spectra and mass functions are shown 
in the top and middle panels of Fig. 9, respectively. We note 
that the cluster mass at the mass function turnover is well- 
preserved, regardless of m„p. The upper limit of the cloud 
mass range controls mostly the high- mass regime of the clus- 
ter mass distribution, that is, the variations with the cluster 
mass of the slope of the mass spectrum, the extent of the 
tail of the mass function above the turnover and the cluster 
mass at the peak of the mass- weighted mass function. 

Burkert & Smith (2000) revisit the issue of how differ- 
ent among galaxies is the high-mass regime of the cluster 
mass distribution. They show that the bright end of the 
luminosity distribution of globular cluster systems in five el- 
hptical and spiral galaxies (M87, NGC5846, NGC4472, the 
Milky Way and M31) is well-fitted with an underlying mass 
spectrum of the form 



dN 
dm 



-3/2 - 

oc m ' e 



m/mc 



(10) 



under the assumption that the cluster mass-to-light ratio is 
constant within each system. In equation 10, the exponential 
cutoff causes the mass spectrum to drop off toward high 
masses more steeply than a pure power-law of the form ^ oc 
^-3/2 rpj^jg approach thus supersedes that of McLaughlin 
& Pudritz (1996), replacing their double index power-law 
description of the mass spectrum high-mass regime (i.e. a > 
—2 and a < —2 below and above rUpeak, respectively) by 
a continuously varying slope. Figure 1 of Burkert & Smith 
(2000) shows that equation 10 fits well the clear curvature of 
the considered globular cluster mass distributions above the 
turnover, the cutoff mass rric being in the range 10® M© < 
mc < 5 X 1O®M0. 

Figure 10 compares the mass distributions of evolved 
cluster systems as obtained by our model (plain symbols) 
to those derived from equation 10 (solid lines, with rric = 
0.5 X 10®, 0.9 X 10®, 2 X 10® Mq in the top panel and rUc = 
1.1 X 10®, 6 X 10® Mq in the bottom panel). For the sake of 
comparison with Burkert & Smith 's (2000) figure 1, we use 
the same cluster mass distribution, namely, the log-normal 
mass function, which is the logarithm of the bell-shaped 
mass function, i.e. log(dA^/dlogm). The bottom and top pan- 
els consider the cases of a power-law mass spectrum and of 
a Gaussian mass function for the protoglobular clouds, re- 
spectively (see section 3.4 for the latter). The evolved cluster 
mass distributions match well the curvature characteristic of 
equation 10, which is the curvature observed for the mass 
distribution of globular cluster systems above the turnover. 
We therefore conclude that our model is able to account for 
the detailed shape of the globular cluster mass distribution 
beyond the turnover. Besides, varying the upper limit m^p 
of the cloud mass range enables us to reproduce different 
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Figure 11. Initial/evolved (solid/open circles) cluster mass dis- 
tributions for three distinct lower cloud mass limits mjo^,. The 
top and bottom panels correspond to the mass spectra and the 
mass functions, respectively. The solid lines are the fiducial case. 
The arrow indicates that the initial mass distributions have been 
shifted by -0.15 in logm, that is, corrected for stellar evolutionary 
mass loss. The evolved mass distributions have been normalized 
so as to contain the same number of clusters as the iiiital ones. 
The initial and evolved distributions are alike, that is, the clus- 
ter mass distribution is preserved over the course of the cluster 
system dynamical evolution. 



exponential cutoff rric at high-masses, while preserving the 
location of the turnover of the mass function, as is observed. 
These variations may originate from a correlation between 
THup and the mass of the parent galaxy and/or its envi- 
ronment. Observations of present-day systems of molecular 
clouds actually reveal a dependence of the upper cloud mass 
limit on the host galaxy (e.g. Wilson et al. 2003, Williams 
& McKee 1997). 



3.3 An equilibrium cluster mass function 

A striking feature of the (practically) universal globular clus- 
ter mass function is its ability to preserve its shape over a 
Hubble-time of dynamical evolution in the tidal field of its 
host galaxy, an effect firstly highlighted by Vesperini (1998) 
(see section 1). The cluster mass distributions predicted by 
our model obey this property, as illustrated in Fig. 11 (top 
panel: mass spectra, bottom panel: mass functions). The 
plain symbols depict the cluster initial mass distributions 
after correction for stellar evolutionary mass loss. The open 
symbols show the corresponding cluster mass distributions 
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Figure 12. Incomplete gamma functions Q of the fit of the 
evolved cluster mass functions to the observed one (thick curves) 
and cluster mass fractions f^Q^^ at an age of 1.3 Gyr as predicted 
by the model (thin curves) versus the mean logmo and standard 
deviation a of the Gaussian cloud mass function. The dotted area 
corresponds to evolved cluster mass functions for which the clus- 
ter mass at the turnover is larger than 3 X 10^ M©. The four- 
branch stars represent different cases presented in Fig. 13 and 
Table 4, the solid symbol being our fiducial case. 



evolved to an age of 13 Gyr and normalized so as to contain 
the same number of clusters as the initial ones. We can thus 
straightforwardly estimate whether the shape of the mass 
distributions evolve significantly over a Hubble-time owing 
to the evaporation and dissolution of globular clusters. The 
initial and evolved cluster mass distributions for the fidu- 
cial case {rriiow = 6 x IO'Mq, solid linos) appear indistin- 
guishable in shape. Slightly different lower cloud mass limits 
{niiow = 2 X lO^M© and miow = lO^M©, dotted and dashed- 
dotted lines, respectively) also result in fairly well-preserved 
cluster mass functions. 



3.4 Another cloud mass distribution 

The bottom panel of Fig. 7 demonstrates that a power- 
law protoglobular cloud mass spectrum with a narrow mass 
range, e.g. lO^M© < rricioud ^ 2 x lO^M©, leads to a good 
fit (Q ~ 0.1) of the modelled cluster mass function onto 
the observed one. This suggests that the present-day halo 
cluster mass distribution, which covers two decades in mass, 
may equally-well arise from a characteristic mass for the 
protoglobular clouds. In order to investigate this point more 
closely, we now assume that the protoglobular cloud mass 
function obeys a Gaussian of mean logmc and of standard 
deviation a. As for the power-law hypothesis presented in 
section 3.1, the outcomes of the simulations are summa- 
rized as the iso-Q and iso-f^Q^^ curves in the (logmc, 
a) plane (sec Fig. 12). We have considered a probability 
distribution p(e) for the star formation efficiency identical 
to that derived in the frame of the power-law hypothesis, 
namely, 6 = —2.9 and rc = 0.025. A Gaussian protoglob- 
ular cloud mass function centered at logmc = 6.15 with a 
standard deviation a less than 0.4 fits well (Q ~ 0.1) the 
Old Halo cluster mass function, while predicting the correct 
halo cluster mass fraction, namely, /"°"'' ~ 2 %. Although 
a higher logmc improves the goodness of fit, it also shifts the 
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Figure 13. Assumed Gaussian cloud mass functions: Cloud mass 
distributions (crosses) and evolved (open circles) cluster mass dis- 
tributions corresponding to some of the cases listed in Table 4. In 
each panel, the solid lines are the fiducial case. The top and mid- 
dle panels show how the evolved cluster mass functions respond 
to variations in the mean logmG a^nd the standard deviation a of 
the Gaussian cloud mass function, respectively. The bottom panel 
displays the evolved mass- weighted mass functions corresponding 
to the cluster mass functions in the middle panel. The Old Halo 
mass distribution is shown as the full squares. 



turnover of the evolved mass function toward unacceptably 
high cluster mass, an efi'ect illustrated by the dotted area at 
the right of the (logmo, cr) plane, for which logmTO > 5.5 
as in Fig. 7. As for the power-law cloud mass spectrum, this 
artifact comes from the slight overestimate of the number of 
low-mass clusters by our model (see section 3.1). 

In Table 4, we investigate how model parameter varia- 
tions affect our results with respect to the adopted fiducial 
case logrriG = 6.15 and a = 0.3 (solid four- branch star in 
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Table 4. Same as Table 3, but for a Gaussian cloud mass function with a mean logmc and a standard deviation cr. 





Innfm /~i 




s 










^ M 


Fk.j 
^ N 


O 


Fiducia.1 case 


6.15 


0.3 


-2.9 


0.025 


0.025 


1.2 X 10~^ 


0.019 


0.28 


0.43 


0.10 


A{log ma) 


5.70 
6.60 


0.3 
0.3 


-2.9 
-2.9 


0.025 
0.025 


0.025 
0.025 


1.2 X IQ-^ 
1.2 X 10"'^ 


0.013 
0.021 


0.19 
0.31 


0.26 
0.53 


10-15 

0.04 


Mfac) 


6.15 
(i.ir, 


0.3 

().:-! 


-1.9 


0.006 

0.018 


0.025 


3.1 X 10"-^ 
()..') X 10 


0.053 
0.008 


0.29 
0.28 


0.46 

0, 10 


0..30 
0.02 


m 


6.15 
6.15 


0.3 
0.3 


-2.4 
-4.4 


0.005 
0.123 


0.010 
0.050 


0.5 X 10"^ 
2.4 X 10"'^ 


0.021 
0.020 


0.29 
0.28 


0.44 
0.41 


0.34 
0.03 


A(a) 


6.15 
6.15 


0.0 
0.6 


-2.9 
-2.9 


0.025 
0.025 


0.025 
0.025 


1.2 X 10"^ 
1.2 X 10-3 


0.018 
0.018 


0.26 
0.27 


0.45 
0.39 


0.13 
0.03 




6.15 


0.3 


-2.9 


0.025 


0.025 


1.2 X 10-3 


0.026 


0.38 


0.61 


0.05 



Fig. 12). The top and middle panels of Fig. 13 illustrate the 
effect of varying logmc and a, respectively. The protoglobu- 
lar cloud mass functions and the evolved cluster mass func- 
tions arc depicted as the crosses and the open circles. The 
shape of the cluster mass function being little affected by dy- 
namical evolution (see section 3.3), the cluster initial mass 
functions are not displayed for the sake of clarity. The main 
controlling parameter of the turnover location is the mean 
protoglobular cloud mass, i.e., the higher logmc, the higher 
logmTO- In case of a lower logmc (e.g. logmc — 5.7), the 
cluster initial mass function favours lower-mass clusters and 
the present-day cluster mass fraction f^^^^ and the sur- 
vival rates Fm and Fn get smaller (see Table 4). As for the 
standard deviation a of the Gaussian cloud mass function, 
it controls the high-mass regime of the cluster mass distri- 
bution, without affecting the turnover location of the cluster 
mass function significantly. A larger standard deviation a re- 
sults in a more extended liigh-niass tail for the cluster mass 
function and a larger cluster mass rUpeak at the peak of the 
mass-weighted mass function (middle and bottom panels of 
Fig. 13, respectively). In that sense, a plays a role similar 
to the upper limit ruup of the power-law cloud mass spec- 
trum (see section 3.2). We have also investigated the effect of 
variations in the star formation efficiency distribution p{e) 
and in the cluster disruption time-scale tdis, while keeping 
the fiducial cloud mass function. The evolved cluster mass 
functions, the survival rates Fm and Fn and the present- 
day cluster mass fractions /q'^"' are very similar to those 
derived in the frame of the power-law hypothesis (compare 
Table 3 and Table 4). 

3.5 The universal globular cluster meiss at the 
turnover and the need for a high 
protoglobular cloud mass-scale 

In order for our model to explain the universal Gaussian 
globular cluster mass function, the protoglobular clouds 
must be characterized by an almost invariant high-mass 
scale, either in the form of a lower truncation of the power- 
law cloud mass spectrum {miow — 6 x lO'^Mg) or in the 
form of a high mean mass logmc ~ 6.2 if the clouds are 
distributed following a Gaussian in log-mass. 

The old age of GCs has prompted the development of 



models for the formation of protoglobular clouds before or 
during the formation of galax;ies, in which the clouds actu- 
ally show a high mass- scale of order 10^ - 10^ Mq . Peebles & 
Dickc (1968) assign a prcgalactic origin to globular clusters 
in that their gaseous precursors were the first gravitationally 
bound systems to form out of the expanding Universe. They 
infer a typical protoglobular cloud mass of order IO^Mq. 
On the other hand. Fall & Rees (1985, see also Kang, Lake 
& Ryu 2000) suggest that globular clusters formed during 
the collapse of the protogalax;y out of gas that cools slowly 
below 10* K, thereby imprinting a preferred mass-scale of 
~ 10® M0 for the protoglobular clouds. Although these char- 
acteristic cloud masses match well those we have derived 
independently in sections 3.1 and 3.4, the prcgalactic and 
protogalactic scenarios of Peebles & Dicke (1968) and Fall 
& Rees (1985), respectively, are hampered with several prob- 
lems. For instance, a prcgalactic origin for globular clusters 
is hardly reconciled with the paucity of intergalactic globular 
clusters, while a metallicity [Fc/H] larger than —2 (i.e. the 
metallicity range of most globular clusters in the Galactic 
halo) prevents the gas temperature from hanging at lO'' K. 
As a result, considerable efforts have gone into analysing 
the protoglobular cloud cooling histories. External sources 
of UV or X-ray flux (a background of hot halo gas, an active 
galactic nucleus or population III stars) have been invoked 
for the gas to retain those characteristic temperature and 
mass (Kang et al. 1990, Murray & Lin 1992). More recently, 
Bromm & Clarke (2002) explored a mechanism for the for- 
mation of the first globular clusters, operating during the 
assembly of dwarf galaxies at high redshift, 2: > 10. Adopt- 
ing the ACDM model for structure formation, they showed 
that the collapse and virialization of a dwarf galaxy result 
in the formation of a handful of high-density gas clumps, 
these being identified as globular cluster progenitors. Be- 
cause, in these models, dwarf galax;ies formed before the 
epoch of reionization (z > 10), no UV background prevents 
tlK> gas temperature from dropping below 10 K and the 
gas thus falls into the dark matter subhalos. The result- 
ing clump masses are determined by the mass spectrum of 
the dark matter substructure. Bromm & Clarke 's (2002) 
three-dimensional numerical simulations of the collapsing 
dark matter and gaseous components lead to a clump mass 
range of 10® < nicioud ^ Mq . Therefore, although the 
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hypothesis that the gaseous progenitors of old globular clus- 
ters are characterized by a mass-scale of order ~ 10® M0, as 
found in our simulations, still remains debated, this appears 
to be a real possibility. 

In contrast, it is clear that prosont-day mass distribu- 
tions of giant molecular clouds and of their cluster form- 
ing cores show no evidence for a low-mass cut-off. Giant 
molecular cloud masses are well described by a featureless 
power- law with a ~ —1.7 down to the completeness limit, 
this being of order 10* Mq in the Large Magellanic Cloud 
(Blitz et al. 2006). If the mass distribution of the cluster 
gaseous progenitors were genuinely truncated at ~ 10'' M0, 
this would load to a boll-shaped cluster initial mass func- 
tion with a cluster mass at the turnover of a few 10'^ Mq 
only. Considering a system of young massive star clusters 
located beyond the Local Group of galaxies whose detec- 
tion limit is at a similar cluster mass, only the high-mass 
regime of the boll-shaped mass function is probed and the 
inferred mass function is a power-law. Therefore, our model 
does not contradict the observations of power-law cluster 
mass spectra, with spectral indices — 2 < a < —1.8 down 
to 10^ Mq — lO'* Mq , for systems of young star clusters in 
present-day starbursts and mergers. ^ 

The exact shape of their mass distribution remains dis- 
puted, however. This is exemplified by the heavily stud- 
ied system of star clusters hosted by the Antennae galax- 
ies (NGC 4038/4039). While Zhang & Fall (1999) report 
a featureless power-law mass spectrum with a slope of —2 
down to 10* Mq, Fritze v. Alvensleben (1998, 1999) finds a 
Gaussian mass function similar to that of old globular clus- 
ters. Because the formation duration of young star cluster 
sytems formed in ongoing or recent starbursts may be a sig- 
nificant fraction of the system's median ago, deriving their 
mass function is hampered by the strong variations in time 
of the integrated cluster mass-to-light ratio. That is, the ob- 
served cluster luminosity distribution may not be a faithfull 
mirror of the underlying cluster mass distribution and reli- 
able cluster age estimates are required to convert the former 
into the latter (Meurer 1995). Moreover, avoiding bright star 
contamination and accounting for completeness effects prop- 
erly prove equally challenging and may lead to discrepant 
cluster luminosity functions from one study to another. For 
instance, Anders et al. (2007) report the first detection of 
a turnover in the Antennae cluster luminosity function at 
Mv — —8.5. The young age of their sample (<100Myr), 
however, implies a cluster mass at the turnover roughly an 
order of magnitude lower than what is observed for old glob- 
ular cluster systems. 

As for the galaxies NGC 3310 and NGC 6745, de Grijs 
et al. (2003) report power-law cluster mass spectra with 
Q ~ — 2. The cluster mass ranges probed corresponding to 
the high mass regime of the globular cluster Gaussian mass 
function (> 10^ Mq and > 4 x 10^ Mq, respectively), one 
cannot distinguish the observed cluster mass functions from 
that of old globular clusters. 

^ That the slope of power-law cluster mass spectra in starbursts 
(—2 ^ a ^ —1.8) is often slightly steeper than that of cluster 
forming cores (a ~ —1.7) -assuming that the value deduced for 
the Galactic disc is valid in other large galaxies- may tell us some- 
thing about an e — mdoud relation. To investigate that particular 
point is beyond the scope of the present study, however. 



Region " B" of the starburst galaxy M82 hosts an al- 
most coeval cluster population. In spite of an intermediate 
age (~1 Gyr), its bell-shaped mass function shows a turnover 
already located at a cluster meiss similar to that of old glob- 
ular cluster systems, i.e. 1.5 x 10'"' Mq. de Grijs, Parmen- 
tier & Lamers (2005) find that that cluster mass function is 
inconsistent with an initial power-law mass spectrum with 
Q ~ — 2, since the initial number of clusters would be un- 
pliysically high in that case ®. In contrast, a Gaussian initial 
cluster mass function similar to that observed today (hence, 
similar to the mass function characteristic of old globular 
clusters) constitutes an acceptable solution. 

An additional concern worth being borne in mind when 
dealing with the mass distribution of young clusters is 
whether they are in virial equilibrium or not. Actually, in 
case of instantaneous gas removal, a protocluster settles 
back into virial equilibrium at an age of 40-50 Myr only (see 
Figs. 2 and 5 of Goodwin & Bastian 2006). That is, any 
population of clusters younger than that age is necessar- 
ily contaminated by unbound clusters whose stellar content 
is being scattered into the field (i.e. those characterized by 
e < eth) as well as by clusters which are still experienc- 
ing star escapes due to gas expulsion (i.e., when e > cth). 
The masses of these out-of-equilibrium/expanding clusters 
are in the ranges - e x m^ioud and Fbound x e x nicioud - 
e X rricioud, respectively. Yet, the mass of relevance here is 
that of gas-free bound star clusters, i.e. Fbound x e x rricioud- 
Consequently, the mass function of clusters younger than 40- 
50 Myr, even if corrected for all the above mentioned caveats 
(temporal cluster mass-to-light ratio variations, bright star 
contamination and completeness effects), is not necessar- 
ily the cluster initial mass function inferred by studies of 
cluster sytcm secular dynamical evolution (e.g., Baumgardt 
1998, Vesperini 1998, Fall & Zhang 2001 and Parmentier 
& Gilmore 2005). For instance, Bik et al. (2003) report a 
power- law cluster mass spectrum with a ~ — 2 and going 
down to 10'' Mq for the cluster population in the inner spi- 
ral arms of the interacting galaxy M51. That mass spec- 
trum is "pre- initial", however, as it covers clusters younger 
than 10 Myr, i.e. clusters still experiencing infant mortal- 
ity and infant weight-loss. Despite being less than 10 Myr 
old, Bik et al. 's (2003) observed cluster mass spectrum 
will stand for the cluster initial mass spectrum which sec- 
ular dynamical evolution models build on, provided that 
both the star formation efficiency e and the star bound 
fraction Ftound arc independent of the cluster gaseous pro- 
genitor mass, as is explicitely assumed in our model. Indi- 
cations that infant mortality actually constitutes a mass- 
independent process, at least for observed cluster masses 
less than 10* Mq, are given in Bastian et al. (2005) for M51 
and FaU et al. (2005) for NGC4038/39 (see also Anders et 
al. 2007). In that featureless power-law mass distribu- 

tion for the gas-embedded clusters is turned into the same 
featureless power-law for the cluster initial mass distribu- 
tion (except for a vertical shift driven by infant mortality 



® This result is driven by the combination of the high cluster mass 
at the turnover and the 1-Gyr age of the clusters. Ongoing HST 
l/-band observations will enable to refine these cluster mass and 
age distributions as well as the initial number of clusters required 
if starting from a power-law. 
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and infant weight- loss ) . Bik ct al. 's (2003) result may thus 
actually points to a featureless power-law initial cluster mass 
spectrum with the "canonical" slope a ~ — 2 in M51. 

The above discussion shows that it is not straightfor- 
ward to derive the initial mass distributions of young star 
cluster systems located beyond the Local Group. In con- 
trast, the Large Magellanic Cloud and the Small Magellanic 
Cloud are close enough so that a detailed survey of even 
faint clusters can be carried out. Building on the framework 
of Boutloukos & Lamers (2003), de Grijs & Anders (2006) 
derive a disruption time-scale log(i4"'/yr) = 9.9 ±0.1 for a 
10* M0 cluster (the disruption time-scale for any other clus- 
ter mass being obtained from t*" = tf^m/10'^MQ)°-'^^, see 
Boutloukos & Lamers (2003) for details). So long a cluster 
disruption time-scale stems from the low-density environ- 
ment of the Clouds and guarantees that the cluster mass dis- 
tributions derived by de Grijs & Anders (2006) have not yet 
boon significantly altered by dynamical secular evolution. 
Considering clusters older than 60 Myr and more massive 
than 10^ M©, the cluster mass spectrum slope is a ~ —2, i.e. 
fully consistent with what is generally inferred for young star 
cluster systems. Similar results are obtained by Hunter et 
al. (2003). Yet, the mass spectrum slope of clusters younger 
than 60 Myr appears significantly shallower, with a ~ —1.7 
for the same mass range and a ~ —1 for clusters less mas- 
sive than 10^ Mq (see Fig. 8 and Table 3 of de Grijs & An- 
ders 2006). They ascribe these substructures in the young 
agc/low-mass regime to ongoing infant mortality. As already 
quoted, at so young an age, the surveyed population of clus- 
ters necessarily consists of a mix of bound and unbound 
clusters. 

The observation of power-law mass spectra for young 
massive star clusters has led to the idea that the initial 
mass spectrum of old globular clusters was also a power-law 
with a similar slope. Yet, when evolving an initial power- 
law, most current models fail to reproduce the proscnt-day 
globular cluster Gaussian mass function, except under very 
specific conditions, generally not consistent with observa- 
tions (e.g. a narrow range of perigalactic distances, see Fall 
& Zhang 2001 vs. Vesperini et al. 2003). In an attempt to 
reconcile the Gaussian mass function of old globular clusters 
with the power-law cluster mass spectrum of young massive 
star clusters, Vesperini & Zepf (2003) build on the observed 
trend between the mass of Galactic globular clusters and 
their concentration (i.e. the more massive the cluster, the 
higher its concentration; van den Bcrgh 1994, see also Fig. 4 
of Larsen 2006). This rough correlation being likely of pri- 
mordial origin (Bellazini et al. 1996), they investigate how 
the dissolution of low-mass low-concentration clusters af- 
fects the temporal evolution of the cluster mass function. 
Their results suggest that it may be possible to reproduce 
the 13 Gyr-old bell-shaped globular cluster mass function 
starting either from an initial power-law mass spectrum or 
from a Gaussian mass function similar to that today. A de- 
tailed study of this effect must still be carried out, however. 

If the protogalactic era actually sets a characteristic 
mass logmo for the protoglobular clouds (or a lower mass 
limit miow in case of a power-law protoglobular cloud mass 
spectrum) independent of the host gala^xy, then our gas re- 
moval model naturally explains why the initial mass function 
of old globular clusters is a Gaussian similar to that observed 



today, while the mass function of present-day massive star 
clusters is a power-law. That the initial mass functions of 
old globular clusters and of present-day massive star clus- 
ters are different is also suggested by the work carried out 
by Hunter ct al. (2003). In their study of about 1,000 star 
clusters of the Large and Small Magellanic Clouds, they de- 
termine the cluster initial mass function by fitting cluster 
population models to the mass distribution integrated over 
age and to the age distribution integrated over mass. They 
derive a power-law cluster mass spectrum with a slope > —2. 
They note however that the distributions of the oldest clus- 
ters, i.e. those corresponding to Galactic halo globular clus- 
ters with respect to their age and mass, do not appear to 
be extrapolations of the distributions for the less massive 
and younger clusters. To explain why the old and massive 
clusters form a distinct population, Hunter et al. (2003) pro- 
pose that they formed under a different cluster initial mass 
function, presumably a Gaussian similar to that of the Old 
Halo cluster system. 

In our model, old globular clusters formed out of ex- 
clusively massive (~ 10® M©) gaseous precursors, span- 
ning a one decade or so mass range (sections 3.1 and 3.4). 
This is significantly narrower than the 4-decade mass range 
of present-day giant molecular clouds and giant molecular 
cloud cores. That narrow mass range is then broadened and 
turned into a bell-shaped cluster initial mass function by 
gas removal, while the ~ 10* Mq protoglobular cloud mass- 
scale guarantees that the turnover settles at the observed 
value. This difference in the mass distribution of cluster 
gaseous progenitors does not necessarily imply markedly dis- 
tinct locations for old globular clusters and young massive 
star clusters in the fundamental plane of star clusters. The 
fundamental plane of star clusters describes how they are 
distributed in terms of mass, compactness and velocity dis- 
persion. Since gas removal leads to the expansion of the pro- 
toclustcr and to the loss of a fraction of its stars, specifically 
those with the hotest kinematics, it certainly plays a key role 
in assigning star clusters to their initial location in the fun- 
damental plane. The high mass-scale of protoglobular clouds 
may result from a different distribution of their sources of 
support against gravitational collapse. That is, the contri- 
butions of thermal motions, turbulence and magnetic fields 
to the protoglobular cloud overall velocity dispersion may 
be different from what is observed today. For instance. Fall 
& Rees (1985) build on the hypothesis that protoglobular 
clouds are thermally supported only, while Harris & Pu- 
dritz (1994) propose that the high mass of the protoglobular 
clouds (the star forming cores of their " supergiant molecu- 
lar clouds") stems from strong non-thermal motions driven 
by magnetic fields and turbulence. While gas removal de- 
pends on gas velocity dispersion, the lower the velocity dis- 
persion, the quicker the gas removal (e.g. Elmegreen & Efre- 
mov 1997, Kroupa & Boily 2002), the source of the velocity 
dispersion itself is of little relevance. As a result, cluster 
forming clouds/cores differing with respect to how they are 
gravitationally supported should not lead to markedly dif- 
ferent initial locations in the star cluster fundamental plane. 
Young massive star clusters do occupy a wider region of the 
so-called K-spacc (i.e. the fundamental plane described in 
terms of cluster mass and cluster compactness) than the 
old globular clusters do, even when the former are aged to 
lOGyr using Simple Stellar Population models in order to 
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account for their stellar evolution-driven fading (see Fig. 8 
of Bastian et al. 2006). This difference is most probably the 
imprint of the cluster dynamical evolution (rather than dif- 
ferences in the formation process), in the sense that only 
a narrow region of the mass-radius plane of globular clus- 
ters is stable against a Hubble-time of dynamical evolution 
(Gnedin & Ostriker 1997). 

While our simulations suggest the existence of a high- 
mass scale for the protoglobular clouds, which the proto- 
galactic era may actually be responsible for, we emphasise 
that our model still lacks a crucial ingredient which pro- 
vents us from drawing any definitive conclusion in this re- 
spect. The Fbound vs. e relation used in our Monte-Carlo 
simulations has been derived for the specific case of isolated 
clusters, that is, it ignores the effect of an external tidal field 
upon the protocluster evolution. Since residual gas expulsion 
leads to protocluster expansion (see e.g. Fig. 7 of Goodwin 
1997, Fig. 1 of Geyer & Burkert 2001), a fraction of oth- 
erwise supposedly bound stars may be driven beyond the 
protocluster tidal radius. That is, an external tidal field un- 
binds a larger fraction of protocluster stars than in the case 
of an isolated group of stars (i.e. the Fbound vs. e relations of 
Fig. 1 define upper limits rather than one-to-one relations). 
At a given galactocentric distance, a low-mass cluster has 
a smaller tidal radius and, therefore, a greater probability 
of additional star loss (with respect to Fig. 1) than a high- 
mass one. The galactic tidal field may thus well lead to addi- 
tional destruction of low-mass clusters. This effect may also 
be enhanced by the fact that low-mass clusters preferen- 
tially arise from low star formation efficiencies, inducing in 
turn a greater spatial expansion of their stellar component 
(see e.g. Fig. 3 of Geyer & Burkert 2001). Detailed A^-body 
simulations are therefore required to investigate whether an 
external tidal field may infiucncc the turnover location. For 
instance, let us consider the case of a small lower limit for 
the cloud mass range, say, logmjo-ui ~ 4 (or a characteristic 
protogobular cloud mass logmc << 10^ M0). In the present 
stage of our model, in which clusters are regarded as isolated 
star systems, the resulting turnover location is predicted at 
logmTO ^ logmjou, ~ 4, more than an order of magnitude 
below what is observed. The question is whether a model ac- 
counting for the galactic tidal field properly may drive that 
low-mass turnover towards higher cluster masses, possibly 
up to the equilibrium location at logmro — 5.3, as a result 
of protocluster tidal truncation. Should this be the case, a 
lower-mass limit miow for the cloud mass range significantly 
less than 6 x lO' M0 (or a characteristic protoglobular cloud 
mass << 10^ Mq), not necessarily universal among galax- 
ies, may actually be acceptable. We note however that if 
an external tidal field contributes significantly to the shape 
of the cluster initial mass function by removing low-mass 
protoclusters and creating a turnover even before secular 
evolution sets in, then the effect should also be observed 
for the mass function of young massive star clusters formed 
in the dense environment of starbursts galaxies. Clearly, we 
need more detailed protocluster gas removal modelling (tak- 
ing into account external effects) as well as more surveys of 
young star cluster systems located beyond the Local Group, 
going deep and correcting for bright star contamination and 
incompleteness effects. 



3.6 The issue of discreteness 

Up to now, each of our Monte-Carlo simulations has encom- 
passed 10^ clouds. With the star formation efficiency distri- 
bution p{e) inferred in section 3.1 {6 = —2.9, rc = 0.025), 
this corresponds to an initial number of clusters of about 
44,000 (i.e. ~ 44,000 clouds achieve e > eth)- In Figs. 7-13, 
the evolved cluster mass functions have been normalised so 
as to contain 72 clusters with 10* < m < 10® M©, namely 
the number of Old Halo clusters in that mass range. The ini- 
tial cluster mass functions have been normalised accordingly. 
With so high an initial number of clusters in each simula- 
tion, discreteness issues inherent to any size-limited sample 
are avoided and the results (cluster mass functions, survival 
rates Fm and Fn and present-day cluster mass fractions 
fcc^^) are independent of the random seeds. Nevertheless, 
Table 3 and Table 4 show that the fraction Fn of surviv- 
ing clusters is < 0.5, implying an initial number of Old Halo 
clusters > 150 only. Therefore, it is interesting to investigate 
how the fiducial cases behave if the number of protoglobular 
clouds is limited so as to lead to a cluster system containing 
as many members as the Old Halo at an age of 13 Gyr. 

This section presents the results of simulations for which 
the number of protoglobular clouds is Ncioud — 43, 000. We 
have reconsidered both fiducial cases, namely, the power-law 
cloud mass spectrum and the Gaussian cloud mass function. 
For each cloud mass distribution. Tables 5 and 6 present 
the goodness of fit Q, the present-day cluster mass frac- 
tion f^Q^^ , the survival rates Fm and Fn and the number 
of clusters Nac at an age of 13 Gyr for a sequence of six 
simulations. The mass range of the Old Halo clusters being 
10'' < m < 10® Mq, the column Nqc quotes both the total 
number of clusters and the number of clusters in the range 
10* < m < 10® Mq. This second number is used to compare 
the final size of the modelled cluster system to that of the 
Old Halo (i.e. 72 clusters) consistently. The corresponding 
cluster evolved mass functions are shown and compared to 
the observed one in Fig. 14. In contrast to the previous sim- 
ulations, no vertical shift has been applied to the modelled 
distributions. Also shown as the thick solid lines arc the 
mass functions derived if Ndoud = 10^ (i.e., no discreteness 
issue, see sections 3.1 and 3.4). 

The successive simulations lead to widely different val- 
ues for the incomplete gamma function, ranging from un- 
acceptable fits {Q = 10~^) up to excellent ones (Q > 0.1). 
Globally however, most fits are acceptable with Q > 0.001. 



3.7 The origin of halo field stars 

In the present-day Galaxy, the formation of unbound stellar 
groups is the rule and not the exception (see section 2.2.2). 
Most field stars in the Galactic disc likely originate from 
embedded clusters which either lost a fraction of their orig- 
inal members or wore disrupted while emerging out of their 
natal clouds. Embedded clusters may therefore well be the 
fundamental units of star formation in the sense that they 
may account for the formation of the vast majority of all 
stars in the Galaxy (Lada & Lada 2003). 

It is likely that that paradigm also characterizes the 
formation of the stellar halo. Actually, 98 % of its mass con- 
sists of field stars and so large a mass fraction cannot be 
accounted for by the secular evaporation and disruption of 
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Table 6. Same as Table 5, but for the fiducial case of the Gaussian 
cloud mass function. 



4.5 5 5.5 
logm/Mg 

Figure 14. Top panel: Evolved cluster mass functions at an age 
of 13 Gyr for the fiducial case of a power-law cloud mass spectrum. 
The thick and thin curves correspond to protoglobular cloud num- 
bers of Ncioud = 10^ a-nd Ndoitd = 43, 000, respectively. In the 
second case, the limited size of the sample results in a different 
mass function for each simulation. The Old Halo mass distribu- 
tion is shown as the full squares. Bottom panel: Same in case of 
an underlying Gaussian cloud mass function. 



Table 5. Effects of sampling noise: Results of different simula- 
tions run with N^ig^^i ~ 43, 000 in the fiducial case of the power- 
law protoglobular cloud mass spectrum. The limited initial num- 
ber of clusters (~ 190) implies that each run provides a different 
cluster mass function, as shown by the thin solid lines in Fig. 14. 
In addition to the values of Q, f^Q^^ , Fm and Fjv, the table also 
lists the number Nqq of clusters which survive over a Hubble- 
time, both total and in the mass range 10* < m < 10® Mq, i.e. 
the mass range of the Old Halo clusters. 



Q 


fac 


Fm 


Fn 


Ngc 


0.068 


0.019 


0.30 


0.46 


(87, 75) 


0.001 


0.016 


0.25 


0.39 


(74, 71) 


0.420 


0.018 


0.29 


0.42 


(79, 70) 


0.010 


0.013 


0.21 


0.38 


(72,68) 


0.226 


0.019 


0.29 


0.43 


(81,71) 


10-5 


0.014 


0.25 


0.45 


(86, 80) 
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Fm 


Fn 
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0.065 


0.019 


0.30 


0.44 


(85, 75) 


2 X 10-"' 


0.018 


0.28 


0.40 


(76, 73) 


0.326 


0.017 


0.27 


0.38 


(72,67) 


4 X 10-5 


0.016 


0.25 


0.47 


(90, 83) 


0.139 


0.019 


0.30 


0.47 


(90, 77) 


0.021 


0.016 


0.25 


0.38 


(73, 68) 



globular clusters over a Hubble-time, regardless of the shape 
of the cluster initial mass function (e.g. Vesperini 1998, 
Baumgaxdt 1998, Parmentier & Gilmore 2005). For instance, 
Tables 5 and 6 demonstrate that the initial total mass in halo 
clusters is at most 10® Mq (see the survival rate Fm), which 
is an order of magnitude lower than the mass of the stellar 
halo. The present-day overwhelming fraction of field stars 
thus constitutes an " ab initio" feature of the Galactic halo. It 
is worth noting that this may help understand why field stars 
and globular cluster stars sometimes show different patterns 
in their light element abundances. For instance, the group of 
O-poor and Na^rich stars observed in some globular clusters 
is hardly detected in the field (Carretta, Gratton & Sneden 
2000). Chemical abundance anomalies observed in globular 
cluster stars are usually ascribed to the accretion onto their 
surface of stellar winds of intermediate mass stars ascend- 
ing the asymptotic giant branch, a process made feasible by 
the dense stellar environment characteristic of globular clus- 
ters. With the possible exception of binary systems, ah initio 
field stars remain unaffected by this external pollution, their 
external layers having never been exposed to any accretion 
process. If most present-day field stars have actually never 
belonged to bound globular clusters, then the absence of 
abundance anomalies characteristic of the asymptotic giant 
branch phase for these stars is not surprising. 

The violent relaxation phase affecting protoclusters fol- 
lowing the expulsion of their residual star forming gas consti- 
tutes a prime candidate to explain the origin of field stars in 
the Galactic halo, without confiicting with the well-accepted 
paradigm following which most stars form in clusters. In our 
fiducial cases, at an age of 13 Gyr, the total mass fraction of 
field stars fps = 98% in the halo arises from the three fol- 
lowing contributions (results for the power-law cloud mass 
spectrum and the Gaussian cloud mass function arc similar): 

(i) the disruption of protoclusters in the e < eth regime 
(infant mortality): fps =91%, 

(ii) the loss of stars by protoclusters for which e > eth 
(infant weight-loss): f%g =4%, 

(iii) the evaporation and the disruption of globular clus- 
ters over a Hubble-time (adult mortality): fpg = 3 %. 

Most halo field stars would thus be given off by star 
forming regions whose efficiency e is less than the star forma- 
tion efficiency threshold tth- In this respect, the large value 
for fpg directly results from the steep slope of p(e), namely, 
5 = —2.9 (combined with the assumption of instantaneous 
gas removal, that is, ttn takes its highest possible value and, 
thus, strongly limits the number of clouds with e > eth)- We 
caution however that this result is utterly model-dependent, 
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having been derived under the assumption of a single func- 
tional form for the star formation efficiency distribution p(e) 
valid over the entire range < e < 1, which may not be true. 



4 DISCUSSION AND CONCLUSIONS 

In this paper, we have presented simulations highlighting 
how the mass function of protoglobular clouds evolves into 
that of gas-free bound star clusters as a result of the ex- 
pulsion of the residual star forming gas due to supernova 
activity. The initial mass of a star cluster depends on the 
mass fraction of its gaseous precursor which is turned into 
bound stars. To retain a bound core of stars, a star forming 
region must achieve a star formation efficiency threshold eth 
of about 33% (although this threshold gets lower if the gas 
is removed non-explosively, see section 2.2.3, and/or if the 
number of stars in the embedded cluster is a few thousands 
only). Above this critical value, the fraction Ftound of bound 
stars steadily increases from (when e ^ eth) up to 1 (when 
e < 1) (see Fig. 1). Therefore, the initial mass of a star clus- 
ter does not depend on the star formation efficiency only. 
It is the bound star formation cfHciciicy Fbound x e which is 
important (equation 1). By virtue of the large variations in 
Abound xe (i.e. between and 1), the cluster initial mass func- 
tion does not necessarily mirror the cluster forming cloud 
mass function. In order to investigate this effect, we have 
convolved cloud mass spectra AN /Am with star formation 
efficiency distributions p{t), each value of the star format 
tion efficiency e being coupled to the corresponding fraction 
Fbound of stars remaining bound following gas removal. We 
have explicitcly assumed that e is independent of the cloud 
meiss (section 2.1). 

Starting the simulations with power-law cloud mass 
spectra truncated at a lower-mass limit miow, we have shown 
that the resulting cluster initial mass functions arc bcU- 
shapcd. Since the gaseous precursors of old globular clusters 
(i.e. those formed during the protogalactic era) may actu- 
ally be characterized by a lower mass limit (section 3.5), 
this result provides support to the early suggestion made 
by Vospcrini (1998) that the globular cluster initial mass 
function is a Gaussian similar to that observed today (i.e. 
an equilibrium cluster mass function, see section 3.3). That 
is, the present-day Gaussian mass function is the preserved 
imprint of the cluster formation process rather than the im- 
print of cluster system dynamical evolution. Consequently, 
our model also provides a new approach to investigate the 
origin of the universal location of the globular cluster mass 
function turnover. 

In this respect, we have discussed how the cluster initial 
mass function responds to variations in the input parameters 
of our model (see section 2.2). We have successively varied 
the slope a of the cloud mass spectrum, its lower and up- 
per limits miovj and m^p, the slope 5 and the scale- length 
rc of the star formation efficiency distribution pit), and the 
efficiency threshold tth required to retain a bound core of 
stars (equivalently the gas removal time-scale Tgr measured 
in units of the protocluster crossing-time Tcross). We find 
that the turnover location is mostly sensitive to the lower 
limit miow of the protoglobular cloud mass spectrum. The 
observed universality of the turnover of the old globular clus- 
ter mass function would therefore originate from a common 



value among galaxies for the lower mass limit of protoglob- 
ular clouds, possibly with second-order variations driven by 
differences in the slopes of the cloud mass spectrum, that of 
the efficiency distribution p(e), as well as by differences in 
the gas removal time-scale. 

Combining our model generating cluster initial mass 
functions with the cluster evolutionary model of Baumgardt 
& Makino (2003), we have subsequently investigated which 
input parameters reproduce both the present-day mass func- 
tion of the Old Halo clusters and the present-day mass frac- 
tion in the stellar halo (section 3.1). Considering the case 
of a power-law mass spectrum for the protoglobular clouds 
with a spectral index a ~ —1.7 (as is often reported for the 
present-day giant molecular clouds and their cores) and ex- 
plosive gas removal (i.e. the Fbound vs. e relation considered 
is the solid curve in Fig. 1), we find a good fit for 5 ~ —2.9, 
rc ~ 0.025, miou, ~ 6 x 10^ M© and m„p > 5 x lO^M©. 
With respect to the cluster mass function, the most criti- 
cal parameter is the cloud mass lower limit miow, since it 
sets the turnover location. The cloud mass upper limit m^p 
affects the detailed shape of the mass function beyond the 
turnover without altering the location of the latter. Varia- 
tions in m„p may thus explain why the high-mass regime 
of the cluster mass function differs among galaxies, even 
though the turnover remains unaffected (see section 3.2). 
Variations in the slope 5 and core Tc of the probability dis- 
tribution p(e) of the star formation efficiency weakly influ- 
ence the shape of the cluster mass function at an ago of 
13Gyr (see e.g. Table 3 and Fig. 8). These parameters are 
mostly used to fix the global star formation efficiency e (i.e. 
the mass fraction of gas turned into stars at the scale of the 
protoglobular cloud system) and the mean bound star for- 
mation efficiency tb (i.e., the mass fraction of gas ending up 
in clustered stars following gas removal) (see the top panel of 
Fig. 7). Specifically, for a given value of e, the slope 5 sets, in 
turn, the fraction of protoglobular clouds achieving e> eth, 
eb and the present-day mass fraction of clusters in the halo 
(see equation 9). The total mass of the stellar halo is dom- 
inated by field stars and this overwhelming mass fraction 
(98 %) cannot be accounted for by the secular dynamical 
evolution of globular clusters. The Galactic stellar halo has 
thus been dominated by field stars since its earliest stages. 
In our model, this feature is accounted for by a steep value 
for the slope of p(e), i.e. 5 = —2.9. This results in the vast 
majority of the protoglobular clusters going into the infant 
mortality regime, that is, they arc dislocated following gas 
removal and their stars are scattered into the field. 

Finally, we have also investigated whether a bell-shaped 
cluster initial mass function similar to that observed in old 
clusters today may arise from a characteristic cloud mass, 
the cluster mass spread resulting from the range in the 
bound star formation efficiency Fbound x e. Actually, an en- 
semble of protoglobular clouds with a mass function centered 
around logmc — 6.1 — 6.2 and a standard deviation a < 0.4 
reproduces the Old Halo cluster mass function, the turnover 
location being mostly driven by logmc (see section 3.4). 

In order to generate a Gaussian globular cluster ini- 
tial mass function similar to that today, our model thus 
requires a protoglobular cloud mass-scale of order 10^ M©. 
The shape of the ~ 1-decade cloud mass distribution (trun- 
cated power-law or Gaussian) is of little relevance, as gas 
removal broadens it and turns it into a bell-shaped cluster 
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initial mass function (although the shape of the cloud mass 
distribution does influence the detailed shape of the cluster 
initial mass function) . The observed location of the turnover 
at a cluster mass of ~ 2 x 10® M© stems from the 10^ Mq 
mass-scale. Therefore, the universality of the turnover of the 
globular cluster mass function would originate mostly from 
a common protoglobular cloud mass-scale among galaxies. 
That such a high mass-scale did actually exist remains de- 
bated (section 3.5). In contrast, present-day giant molecular 
clouds and their cores show no such cut-off at high mass. If 
their power-law mass spectrum is truncated at, say, 100 M©, 
the turnover generated by our model gets located at similar 
a cluster mass and the inferred cluster initial mass spectrum 
is a power-law down to < 100 Mq . 

Therefore, with the mass-scale of the cluster gaseous 
progenitors as its key-parameter, our model can account for 
both the Gaussian globular cluster initial mass function and 
the observed power-law mass spectrum of young massive 
star clusters formed in starbursts and mergers. 

We emphasize however that the gas-removal phase of 
our model does not include the impact of external effects, 
e.g., the tidal field of the host galaxy. 

While the mass function of globular clusters is a bell- 
shape, that of open clusters is a power-law (van den Bergh & 
LaFontaine 1984, Battinelli, Brandimaxti, Capuzzo-Dolcetta 
1994). As for massive young clusters formed beyond the 
Local Group, the difference may be driven by the pres- 
ence or not of a lower mass limit for the cluster forming 
clouds/cores. Additionally, we caution that the application 
of our model is necessarily restricted to the initial mass func- 
tion of massive star clusters. Actually, as for low-mass clus- 
ters (e.g. open clusters), their relaxation time-scale is signifi- 
cantly shorter by virtue of their smaller rmmbcr of stars. As 
a result, the gravitational interactions between stars con- 
tribute to cluster dynamical evolution at an earlier time, 
while the exposed cluster is attempting to settle to a new 
state of equilibrium. This eases the formation of a bound 
core of stars, thereby altering the Fbound vs. e relation with 
respect to what is shown in Fig. 1. For instance, despite 
applying a tidal field and despite considering instantaneous 
gas removal, Kroupa, Aarseth & Hurley (2001) find that a 
substantial bound core of stars Ftound — 0-33 remains even 
if e is as low as ~ 0.33. Accounting for the variations in the 
bound fraction of stars as a function of the protocluster mass 
will eventually lead to a picture of the overall cluster mass 
function, from the low-mass regime (100 - 1,000Mq) up to 
the high-mass one (IO'^-IO^Mq) (Kroupa & Boily 2002). 

This paper has been typeset from a T^JC/ I^T^ file prepared 
by the author. 
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